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Abstract: 

^ . We calculate the main observables in B^^^ {p,u,K*)'y and Bg — > (_ft'*,0)7 decays, 
i.e. branching ratios and CP and isospin asymmetries. We include QCD factorisation 
results and also the dominant contributions beyond QCD factorisation, namely long- 
distance photon emission and soft-gluon emission from quark loops. All contributions 
beyond QCD factorisation are estimated from light-cone sum rules. We devise in partic- 
ular a method for calculating soft-gluon emission, building on earlier ideas developed for 
analogous contributions in non-leptonic decays. Our results are relevant for new-physics 
searches at the B factories, the LHC and a future super-flavour factory. Using current 
experimental data, we also extract |Vfd/Vts| and the angle 7 of the unitarity triangle. We 
give detailed tables of theoretical uncertainties of the relevant quantities which facilitates 
future determinations of these CKM parameters from updated experimental results. 
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1 Introduction 



The flavour-changing neutral-current (FCNC) transitions 6 — > 57 and b ^ d'j are among 
the most valuable probes of flavour physics. Assuming the Standard Model (SM) to be 
valid, these processes offer the possibihty to extract the CKM matrix elements |T4(d,s)|, in 
complementarity to, on the one hand, the determination from B mixing and, on the other 
hand, the SM unitarity triangle (UT) based on the tree-level observables and the 

angle 7. These decays are also characterized by their high sensitivity to new physics (NP) 
contributions and by the particularly large impact of short-distance QCD corrections, see 
Ref. [1] for a review. Considerable time and effort has gone into the calculation of these 
corrections which are now approaching next-to-next-to- leading order accuracy [2, 3]. On 
the experimental side, both exclusive and inclusive 6 — > 57 branching ratios are known 
with good accuracy, 5% for B — > K*^ and 7% for B — > X57, while the situation is less 
favourable for b ^ d'y transitions: measurements are only available for exclusive channels. 
In Tab. 1 we give the branching ratios of all established exclusive b — > {d, 5)7 channels. 

Whereas the inclusive modes can be computed perturbatively, using fixed-order heavy 
quark expansion or soft-collinear effective theory (SCET) [2, 7], the treatment of exclu- 
sive channels is more complicated. With presently available methods it is impossible to 
simulate the full amplitude on the lattice, the reason being the occurence of non-local 
correlation functions associated with the insertion of the electromagnetic interaction op- 
erator into the effective Hamiltonian for 6 — > {s,d)^. Instead, one has to resort to effective 
field theory methods, which yield an expansion in inverse powers of the b quark mass, mi,. 
It was shown, in SCET, that the relevant hadronic matrix elements factorise to all orders 
in as and leading order in l/m^ and can be written as [8] 



{V^mB) = e* 



x\l + 0{l/mb] 



Jo . , , 

;i) 

This formula coincides with that obtained earlier in QCD factorisation (QCDF) to next-to- 
leading order in as [9, 10, 11, 12, 13, 14]. In (1), is the photon's polarisation four- vector, 
Qi is one of the operators in the effective Hamiltonian for b {s, d) transitions, 
is a 5 — > F transition form factor, and 0^, (jj^.y are leading- twist light-cone distribution 
amphtudes (DAs) of the B meson and the vector meson V, respectively. These quantities 
are universal non-perturbative objects and describe the long-distance dynamics of matrix 
elements, which is factorised from the perturbative short-distance interactions included 
in the hard-scattering kernels T/ and T^^ . B K7 decays have also been investigated in 
the alternative approach of perturbative QCD factorisation (pQCD) [15]. 

Eq. (1) is sufficient to calculate observables that are dominated by the leading-order 
term in the heavy-quark expansion, like B{B K*'-/). For B{B {p,uj)'~f), however, 
power-suppressed corrections play an important role, for instance weak annihilation (WA) 
which is mediated by a tree-level diagram. In this case, the parametric suppression by one 
power oi l/rrib is alleviated by an enhancement factor 27r^ relative to the loop-suppressed 
contributions at leading order in l/rrib. Power-suppressed contributions also determine 
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B X 10^ 


BaBar [4] 


Belle [5] 


B X 10^ 


HFAG [6] 


B ip,uj)j 
B+ p+7 

B° ^ p% 

B^ ^u-f 


1.251°:^^ ± 0.09 

i.io^°:i±o.o9 

0.79lgJ^ ± 0.06 
< 0.78 


1 00+0.34+0.10 
-'-■"^^-0.31-0.09 

n p;p;+0.42+0.09 
'J-0. 36-0.08 

1 OCT +0.37+0.07 
-'-•^'-'-0.33-0.06 

n nfi+0.34+0.05 

u.yo_o.27-o.io 


B+ K*+-f 
50 ^ ^*o^ 


40.3 ±2.6 
40.1 ±2.0 



Table 1: Experimental branehing ratios of exclusive b (d, 5)7 transitions. All entries are CP 
averaged. The first error is statistical, the second systematic. B — > {p,Lo)'y is the CP average of 
the isospin average over p and uj channels: 

B{B ^ {p, a;)7) = i ^ ^±7) + ^ \B{B^ ^ A) ± B{B^ - ^7)] } • 

the time-dependent CP asymmetry in S — > I/7, see Refs. [16, 17, 18, 19], as well as 
isospin asymmetries [20] — all observables with a potentially large contribution from NP. 
The purpose of this paper is to calculate the dominant power-suppressed contributions to 
(1) and the resulting branching ratios and CP and isospin asymmetries, using the most 
up-to-date hadronic input parameters for form factors and light-cone DAs. Although 
l/rrih effects arc, in principle, accessible in SCET, the vast majority of studies in this 
framework only includes leading-order effects, the reason being a proliferation of new 
effective operators at power-suppressed accuracy, sec Rcf. [21], whose matrix elements 
induce subleading form factors and are largely unknown. For this reason, in this paper we 
adopt a different approach not based on SCET and calculate power-suppressed corrections 
using the method of QCD sum rules on the light-cone (LCSRs). The present paper is 
an extension of our previous work Ref. [22], where we calculated the ratio of branching 
ratios B{B {p,uj)^)/B{B K*'^) in order to determine the ratio of CKM matrix 
elements |14d/Vts| from data. Some of these power corrections, namely those related to 
WA contributions and the isospin asymmetry in B ^ K*'~f, have already been calculated 
in QCDF [12, 13, 14, 20]. Other power corrections cannot be calculated in the framework 
of QCDF. The most relevant of these come from soft-gluon emission from quark loops 
and long-distance photon emission from soft quarks. We have already calculated some 
of these contributions before, using LCSRs: long-distance photon emission in Ref. [23] 
and soft-gluon emission from charm loops in Ref. [19], using heavy-quark expansion in 
powers of l/rric- In this paper, we complete these calculations and develop a method 
to also calculate soft-gluon emission from light-quark loops, thus allowing us to predict 
branching ratios and isospin and CP asymmetries for exclusive B V'j transitions with 
increased precision. We also include the S,, decays Bs (p'y, which is a 6 — S7 transition, 
and Bg K*'~f, which is 6 — > d'-f. All these decays will be studied in detail at the LHC, 
and those of Bu^ at future super- flavour factories [24]. 

Our paper is organized as follows: in Sec. 2 we introduce notations and recall QCDF 
formulae. In Sec. 3 we calculate the WA contributions and in Sec. 4 the long-distance 
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contributions to the B — » V'j amplitude from heavy- and hght-quark loops. In Sec. 5 
we present results for branching ratios and asymmetries; we summarize and conclude 
in Sec. 6. The appendix contains a discussion of the longitudinal and transverse decay 
constants of vector mesons. 



2 Framework and Basic Formulae 

The effective Hamiltonian for b — > transitions, with D — s,d, reads: 



u 



U=u,c 



CiQY + C2Q^ + CiQi 



1=3.. .8 



(2) 



where A^; = V^jj^Vub- This Hamiltonian implicitly relies on the SM unitarity relation, 
also referred to as the GIM mechanism, 

Ar + Ap) + Ai^) = 0, (3) 

which enters the calculation of the penguin contributions. The operators are given by 

= iDiUj)v-AiUjbi)v-A , = {DU)v-A{Uh)v-A , 

Qs = {Db)v-A ^{qq)v-A , Qa = (A^'j)y-A ^{qjqi)v-A , 

Qb = {Db)v-A ^{qq)v+A , Qe = (A&j)y-A ^iqjqi)v+A , 

9 9 

Q7 = Da^'^il + 75)F^, b + ^tjid Da^\l - lb)F^. b , 

= ^m,Da^%l + ^,)G^,b + ^mDDa^'''{l-^,)G^,b, (4) 

where {qQ)v-A{^R)v±A = (?7/^(l — 75)<5)(^7'*(1 i Ts)-^)- The sign conventions for the 
electromagnetic and strong couplings correspond to the covariant derivative — d^-\- 
ieQfA^ + igT"^Al^. With these definitions the coefficients Cy^g are negative in the SM, 
which is the choice generally adopted in the literature. The above operator basis, which 
we shall label BBL after the authors of Ref. [25], is the same as that of Ref. [12], except that 
Qi and Q2 are exchanged. For some applications, in particular calculations of inclusive 
b — > Dj transitions, a different operator basis proves more suitable: the basis adopted for 
instance in Refs. [26, 27], labelled CMM in the following, has Q^^-^^ = Q^^f, but differs 
in Qi...6- It turns out that we need the Wilson coefficients in both bases: 

• Cf^Q{fi ~ mb), calculated according to Ref. [25], for power-suppressed (WA and 
soft-gluon emission) contributions; 
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L/2 [rrib) 


L/g [mh) 






Gg [rrih) 




— n '^99 


1 nriQ 




— n 874 




— n nni 


— n 3nQ 


\Jllb) 


L 2 [nib) 


c. [nil,) 


r'I3BLY,,,,A 

(---i [lllb) 






'•-8 


-0.189 


1.081 


0.014 


-0.036 


0.009 


-0.042 


-0.170 
















-0.288 


1.133 


0.021 


-0.051 


0.010 


-0.065 


-0.191 



Table 2: NLO Wilson coefficients to be used in this paper, at the scales rrib = 4.2 GeV and 
lih = 2.2 GeV. The coefficients labelled BBL correspond to the operator basis of Ref. [25] and 
given in Eq. (4), whereas CMM denotes the basis of Ref. [26]. We use as{mz) = 0.1176 [28] and 
mt{mt) = 163.6 GeV [29]. Note that Cf^^ and C^^^ are exchanged with respect to the basis of 
Ref. [12] and that Cf^^ = see text. Following [13], the CMM set is used for calculating 

hard-vertex corrections to the QCDF formulae and the BBL set at the lower scale is used 
to calculate hard-spectator corrections. The BBL set at scale is used for the calculation of 
power corrections. 

• ^l!^^(a* ^ f^b): calculated according to Ref. [26], for hard- vertex corrections in 
QCDF which are given in terms of two-loop matrix elements for b — > transitions 
obtained in Ref. [27], in the CMM basis; 

• Cf^kil^ ~ 2GeV) and C^^^^{iJ - 2 GeV) for hard-spectator corrections in QCDF; 
although these coefficients refer to a different basis, it is correct to use them together 
as the corresponding operators Q^^s are identical in both bases and independent of 
the basis chosen for the four-quark operators. 

Numerical values of all Ci are given in Tab. 2. Note that the question whether to use 
Wilson coefficients (and other scale-dependent hadronic quantities) at LO or NLO ac- 
curacy is actually non-trivial. Strictly speaking, NLO accuracy is mandatory only for 
C7, as the hadronic matrix element for this term only is also known to NLO accuracy, 
see below. We will evaluate all 0{as) and power-suppressed corrections using both LO 
and NLO scaling for Wilson coefficients and hadronic matrix elements and include the 
resulting discrepancies in the theoretical uncertainty. 

The calculation of the decay amplitudes of exclusive B V'-f decays also requires 
the knowledge of hadronic matrix elements of type {V^\Qi\B). A complete calculation of 
these quantities is not possible to date, but the leading term in an expansion in l/rrih is 
obtained from QCDF, see Eq. (1). The factorisation formula is valid in the heavy-quark 
limit — > 00 and is subject to corrections of order Aqcd/^6- Some of these corrections 
are numerically very relevant: for instance, the contributions from all operators but Q2 
are loop suppressed; hence, the tree-level WA diagram in Fig. 1, which is suppressed by 
one power of m^, comes with a relative enhancement factor ~ Att^. This contribution. 
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(a) 



(b) 



Figure 1: (a): WA diagram. The square denotes insertion of the operator Q,. Photon emission 
from hues other than the B spectator is power-suppressed, except for emission from the final- 
state quark lines for the operators Q^fi, denoted by crosses, (b): soft-gluon emission from a 
quark loop. Again the square dot denotes the insertion of the operator Qi. There is also a 
second diagram where the soft gluon is picked up by the B meson. 

with the operator Q2, is doubly Cabibbo-suppressed for b ^ S'y transitions, but carries no 
CKM suppression factor in b d'-f transitions. WA can also be induced by the penguin 
operators Q3,„q and in this case carries no CKM suppression in b —>■ 37, but comes 
with small (loop-suppressed) Wilson-coefficients. Other examples for relevant power- 
suppressed corrections are CP and isospin asymmetries, Refs. [16, 19, 20], which actually 
vanish in the heavy quark limit. This indicates that in B —>■ V'j transitions simple 1/mb 
counting is, in general, not sufficient to determine the numerical relevance of a particular 
contribution, but that all relevant factors, 

• order of power-suppression in l/rrih; 

• loop suppression or tree enhancement; 

• CKM suppression; 

• size of hadronic matrix elements; 

have to be taken into account. This is a consequence of the fact that in radiative transitions 
the "naively" leading term in Qj is loop suppressed, which is qualitatively different from 
other applications of QCDF, for instance in tt^tt'^, where the leading hadronic 

matrix element describes a tree-level process. 

The exclusive B process is actually described by two physical amplitudes, one 

for each polarisation of the photon: 

^i(fi) = ^(5 ^ K7i(^)) , ^L(ij) =^(5^y7i(^)), (5) 
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where B denotes a [bq) and V a {Dq) bound state. ^ In the notation introduced in Ref. [12] 
in the context of QCDF, the decay amphtudes can be written as 



V2 



Mr) = ^ {iKra^7RiL)(V) + (Xfra^,niL)(V)) {V^^n)\{Q?'^m . (6) 

The Oy" calculated in Refs. [12, 13], coincide, to leading order in l/rrib, with our a^^, 
whereas set zero in [12, 13]. Our expression (6) is purely formal and does not 

imply that the a^^(^) factorise at order l/m^. As a matter of fact, they don't. The 
operators Q^^^^ are given by 

and generate left- (right-) handed photons in the decay b — > D7. The matrix elements in 
(6) can be expressed in terms of the form factor Tf^^ as 

{V{p,vhLiR){q,e)\Q^;^''^\B) 

= -^m,rf-^(0) [e^'''''^elr):p,q,±i{{e*r)*)ipq) - {e*p)iv*q)}] 



27r2 



mbTf^''{0)SLiR), 



{Vip,rjh,iR)iq,e)\{Q?''^m = --^ m,rf-^(0)5^«) , (7) 

where Sl^r are the hclicity amplitudes corresponding to left- and right-handed photons, 
respectively, and e^^r]^) is the polarisation four- vector of the photon (vector meson). The 
definition of can be found in Ref. [30]; our convention for the epsilon tensor follows 

that of Bjorken & Drell: Tr[7"7'^7"'7^75] = Aie"'^'^^. Up-to-date values for all decays 
studied in this paper are given in Tab. 3. The non-perturbative parameters are taken 
from experiment, where available, from lattice (/b), from QCD sum rules (a^) and from 
QCD sum rules on the light-cone (Ti); for the decay constants /-*- results are available 
from both lattice and QCD sum rules; they are discussed in the appendix. No lattice 
results are available for and only partial results for Ti [44]. The numbers in Tab. 3 



-'^Notc that in this paper K* is a (sq) bound state, in contrast to the standard labelhng, according 
to which K*^ = (ds) and = [sd). This is because the calculation of form factors and other matrix 
elements involves light-cone DAs of the vector meson V and in the standard notation used in that context, 
K* always contains an s quark, and K* an s quark. This distinction is relevant because of a sign change 
of G-odd matrix elements under {sq) <-> (gs), see Tabs. 3, 5, 6. 
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differ slightly from those given in Ref. [30] because we include updates of the hadronic 
input parameters. We do not include isospin breaking in the form factors since it is caused 
by the difference of quark masses and electric charges and expected to be of the order 
of 1% only. This is indeed the size of isospin-breaking in the form factor indicated by 
recent measurements of {K-,n-)e+u and D+ {K^,n^)e+u at CLEO [45]. At 

this point we would also hke to comment on the UT angle 7. The value given in Tab. 3 
comes from Belle's Dalitz-plot analysis of the CP asymmetry in B~ — > (X^tt+tt")!^", 
with K^n+n- being a three-body final state common to both and D°. This method 
to measure 7 from a new-physics free tree-level process was suggested in Ref. [46] and 
has been implemented by both BaBar [47] and Belle [32], but the BaBar result currently 
suffers from huge errors. Other determinations of 7 from QCDF or SCET analyses, or 
SU(3) or U-spin fits of non-leptonic B decays, or global UT fits, all come with theoretical 
uncertainties and/or possible contamination by unresolved new physics, so we decide to 
stick, as a reference point, to the tree-level result of Belle. For all observables with a 
pronounced dependence on 7, i.e. 6 — > ^7 branching ratios and isospin asymmetries, we 
will present results as a function of 7. 

As for the oy coefficients, the o^lO^) QCDF, of order 1 in a 1/m;, expansion, 

a,t{V)=Cr + 0{as,l/mb), (8) 

whereas dj'jiiV) are of order l/rrib [17, 18]. It proves convenient to split these coefficients 
into three contributions which we will investigate separately: 

c^UV) = av^i?'^^'!^) + a^nV) + a^^^y) + . . . , (9) 

where a^^*"°^ is the leading term in the l/rrib expansion; all other terms are suppressed 
by at least one power of m;,. We only include those power-suppressed terms that are 
either numerically large or relevant for isospin and CP asymmetries. The dots denote 
terms of higher order in as and further l/rrib corrections to QCDF, most of which are 
uncalculable. Explicit formulae for a^f^'~^^^, complete to 0{as), can be found in Ref. [13]. 
^f/^ann gjj(.Q(jgg ^j^g 0(1/ nib) Contribution of the WA diagram of Fig. 1 which drives the 
isospin asymmetries and has been calculated, in QCDF, and to leading order in ag, in 
Refs. [13, 20] for p and K*, and in Ref. [14] for lo. Preliminary results for the 0{as) 
corrections to WA in S — > p7 were presented in Ref. [48]. This contribution is also 
relevant for the branching ratio of B ^ {p,uj)^] in this case, also long-distance photon 
emission from the soft B spectator quark, which is 0(l/m^), becomes relevant and has 
been calculated in Refs. [23, 49]. We discuss the WA contributions in Sec. 3. The last 
terms in (9), cij^j^y encode soft-gluon emission from a (light or heavy quark) loop as shown 
in Fig. 1 and are particularly relevant for the CP asymmetry; they will be discussed in 
Sec. 4. In Ref. [20] also another class of l/rrib corrections to S — > K*^ was calculated, 
namely 0{as) corrections to the isospin asymmetry in this decay. As these corrections 
break factorisation (require an infra-red cut-off in the momentum distribution of the 
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CKM parameters and couplings 



A [28] 


\Vc,\ [31] 


\Vub\ 


7 [32] 


asimz) [28] 


q; 


0.227(1) 


42.0(7) X 10-^ 


4.0(7) X 10-=^ 


(53 ±20)° 


0.1176(20) 


1/137 


B parameters 


fB, [33] 


fB. [33] 


^bM) [22] 


AB,()Ufe) 


fJ'h 




200(25) MeV 


240(30) MeV 


0.51(12) GeV 


0.6(2) GeV 


2.2 GeV 




p parameters 


fp 


J p 


ai{p) 


4 (p) 


Tf^^(O) 




216(3) MeV 


165(9) MeV 





0.15(7) 


0.27(4) 




cu parameters 


f. 


J u> 




4H 






187(5) MeV 


151(9) MeV 





0.15(7) 


0.25(4) 




K* parameters 


fK* 


fx* 


ai{K*) [34] 




Tf'^'^*(0) 


Tf^^^*(0) 


220(5) MeV 


185(10) MeV 


0.04(3) 


0.15(10) 


0.31(4) 


0.29(1) 


(f) parameters 


f4> 




4(0) 


4(0) 


Tf^^'^(O) 




215(5) MeV 


186(9) MeV 





0.2(2) 


0.31(4) 




quark masses 


m,(2GeV) [35] 


■mb{mb) [31] 


mc{mc) [36] 


mt(mt) [29] 


100(20) MeV 


4.20(4) GeV 


1.30(2) GeV 


163.6(2.0) GeV 



Table 3: Summary of input parameters. The value of \Vub\ is our own average over inclusive and 
exclusive determinations and the result from UTangles, see Refs. [6, 37, 38]. None of our results 
is very sensitive to \Vui,\- For an explanation of our choice of the value of the UT angle 7, see 
text. As is the first inverse moment of the B meson's light-cone DA. A^^ is obtained from A^, as 
described in the text. The vector meson decay constants /y, fy are discussed in the appendix; 
the values of the Gegenbauer moments aj- are compiled from various sources [22, 39, 40, 41] and 
include only small SU(3) breaking, in line with the findings for pseudoscalar mesons [42]. The 
form factors Ti are obtained from LCSRs and are updates of our previous results [30] , including 
the updated values of the decay constants /p,a;,(/) and of ai{K*) [34, 43]. Note that ai{K*) refers 
to a {sq) bound state; for a (qs) state it changes sign. All scale-dependent quantities are given 
at the scale p = 1 GeV unless stated otherwise. 
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valence quarks in the K* meson) and are numerically small, wc do not include them 
in our analysis. As for a^^, the dominant contributions to a':j^{K*) were calculated in 
Ref. [19]. Here we extend this analysis to other vector mesons and also develop a method 
to include light-quark loops. 

We conclude this section by providing explicit results for the QCDF contributions 
to a/, using the formulae of Ref. [13] and the Wilson coefficients and hadronic input 
parameters collected in Tabs. 2 and 3. The decay constants fv and fy are defined as 

{0\q'y^D\V{p)) = mvfvV,, {0\qa,,D\V{p)) = t{v,Pu - P,Vu)fv ; (10) 

their numerical values are discussed in the appendix. The other parameters in Tab. 3 
pertaining to vector mesons are and 02 , the first and second Gegenbauer moments of 
their transversal light-cone DAs of leading twist. In this paper we do not want to go into 
any detail about light-cone DAs, their conformal expansion in Gegenbauer polynomials 
and the dependence of the Gegenbauer moments on the renormalisation scale, but simply 
refer to the relevant literature [39, 40, 41]. 

It turns out that, at the level of two significant digits, all a^^^^^ are equal and so are 
the a"]^^°^. For central values of the input parameters of Tab. 3 we obtain 

^c,QCDF^^^ = -(0.41 + 0.03i) - (0.01 + O.Oli) , 

a^'^^^^{V) = -(0.45 + 0.07i) + (0.02 - Oi) . (11) 

Here we have split the result into contributions from vertex corrections (first term) and 
hard-spectator interactions (second term). The size of the hard-spectator corrections is 
set by the factor 

9 mBTf^^(0)AB 

Note that our value of Xb^, the first inverse moment of the twist-2 i?-meson light-cone 
DA, is obtained from that of A^^ by a simple scaling argument: 

— (Aqcd + = - — Aqcd , 

which follows from the assumption that the Bq DA peaks at the spectator momentum 
k+ = Aqcd, whereas that of Bg peaks at Aqcd + ^s- 

The parameters a^'^'^^^, at tree level, were obtained in Ref. [19] and read: 

0'7R {K A) = C7 — , ay'£ {p,u;,K )^C7 — . (13) 

3 Weak Annihilation Contributions 

The intrinsic WA diagram is shown in Fig. 1; the weak interaction operator is one of the 
charged-current or QCD-penguin operators. All these contributions are 0(l/mb); photon 



hv-^-^ ^^ S-il... ■ (12) 
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WA 


B- K*- 


^0 _^ J^*0 


B^{p,u;) 






induced by 


C (and P) 


P 


C and P 


P 


P 


CKM 


(and 1) 


1 


1 


1 


1 



Table 4: Parametric size of WA contributions to B ^ Vj. C denotes the charged-current 

operators Qi,2) P the penguin operators <53,4,5,6j their Wilson coefficients are small, see Tab. 2. 
CKM denotes the order in the Wolfenstein parameter A with respect to the dominant amplitude 
induced by Qr. 



emission from the b quark and the quarks in the vector meson is further suppressed and 
0(l/m^) - unless the weak interaction operator is Qq^, which can be Fierz transformed 
into {D{1 + 75)g)(^(l — 75) and picks up an additional factor from the projection 
onto the B meson DA which results in this contribution being 0{l/mi,). In Tab. 4 
we show the relative weights of these diagrams in terms of CKM factors and Wilson 
coefficients. The numerically largest contribution occurs for B^ — * P^7- it comes with 
the large combination of Wilson coefficients C2 + C1/3 = 1.02 and is not CKM suppressed. 
For 5° — i> (p°,a;)7 it comes with the factor Ci + C2/3 = 0.17 instead and an additional 
suppression factor 1/2 from the electric charge of the spectator quark {d instead of u). For 
all other decays, WA is suppressed by small (penguin) Wilson coefficients. We evaluate 
the annihilation diagrams at the scale /i — nib. Apart from B — > (p, a;)7, WA is not 
relevant so much for the total values of q^l, but rather for isospin breaking, which is set 
by photon emission from the spectator quark. WA is the only mechanism to contribute 
to isospin asymmetries at tree- level; see Ref. [20] for 0(0;^) contributions. 

Formulae for a^£^'^"(p, K*) in QCDF can be found in Refs. [13, 14]; in this approxima- 
tion, there is no contribution to 07^^'"'. In QCDF, the a^^"^"' are expressed in terms of the 
hadronic quantities 

,y 2%^ fBmvfv w 47r2 fsf^ f\ 0^y(^^) 



" Tf-^(O) mBrribXB ' ^ Tf-^(O) m^m, I v ^'^> 

and dy , obtained by replacing 1/v 1/v in the integrand; is the twist-2 DA of 

a transversely polarised vector meson. For Bg decays one has to set fs fss S'-nd 
correspondingly for the other B meson parameters. Numerically, one finds, for instance 
for the p, hf = 0.22 and dP = -0.59, at the scale p = 4.2 GeV. As Ti ~ l/mf^ and 
/b ~ ^h^^'^ in the heavy quark limit, these terms are 0(l/mb), but not numerically small 
because of the tree-enhancement factors of tt^. 
For uj, K* and (f) we obtain 



a 



u,ann. 



^)Iqcdf = Qdb'^iai + 2(03 + as) + 04) + + d^)ae , 



7L V"^; I QCDF 



«7r°(^)lQCDF = Qdb''('2(a3 + as) + 04) + g<i« + d'^)ae , 
aSr""(0) = Qsb^as + as) + Qs{dt + dt)a^ , 

QCDF 
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(a) (b) (c) 

Figure 2: Example radiative corrections to weak annihilation. The corrections to the B vertex 
in (a) are known [50] and those to the V vertex in (b) are included in /y. For the non-factorisable 
corrections in (c) only preliminary results are available, see text. 



a^r"(^*) = Qsh"" a, + Q,{d^ Q,IQ, + )a, , (15) 

QCDF 

with ai = Ci + C2/3, 03 = C3 + C4/3, 04 = C4 + C3/3, as = C5 + Ce/S, ag = Cq + Cg/S; 
note that ai ^ 02 as compared to [13] as in our operator basis (i.e. the BBL basis) Qi 
and Q2 are exchanged. The expressions for and K* are new; for u, we do not agree 
with [14]. For completeness, we also give the annihilation coefficients for the p, K* and 
a;, as obtained in Ref. [13]: 

«7r^(i^*°)|QCDF = + a,{dr + dT)] , 

«7r°(P°) Iqcdf = Qd + «46'' + a^{dP + dg)] , 



a7r"(P") Iqcdf = Q« [a^^'" + + a^{Qd/Qud', + dg)] , 



«7r(^*°) Iqcdf = ["46^* + a,{dT + *)] , 
4'r"(^*") Iqcdf = Qn [046^* + a6(g./Q„df + rff *)] , 
«7r(P°) Iqcdf = Q<i \.^^^' + ^eC^S + 4)] , 

«7f'(P") Iqcdf = + «6(gd/g.< + <)] • (16) 

Apart from p and a;, all these coefficients are numerically small and do not change the 
branching ratio significantly; the terms in oe, however, are relevant for the isospin asym- 
metries. 

In view of the large size of a^'l™^{p) it is appropriate to have a look at further cor- 
rections. The most obvious ones are 0{as) corrections to the QCDF expressions, shown 
in Fig. 2. As it turns out, the corrections to the B vertex in Fig. 2(a) are known: they 
also enter the decay B — > 7^1/ and were calculated in Ref. [50]. Numerically, they are 
at the level of 10%. Fig. 2(b) shows the vertex corrections to the V vertex, which are 
actually included in the decay constant /y. For the non-factorisable corrections shown 
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in Fig. 2(c) preliminary results have been reported in Ref. [48]; according to [48], these 
corrections are of a size similar to the B vertex corrections. Another class of corrections 
is suppressed by one power of rrib with respect to the QCDF contributions and is due 
to long-distance photon emission from the soft B spectator quark. A first calculation of 
this effect was attempted in Ref. [49] and was corrected and extended in Ref. [23]. The 
long-distance photon emission from a soft-quark line requires the inclusion of higher-twist 
terms in the expansion of the quark propagator in a photon background field, beyond the 
leading-twist (perturbative) contribution; a comprehensive discussion of this topic can be 
found in Ref. [51]. The quantity calculated in Ref. [23] is 

(P~ (p)7(?) I (du)v-A(ub)v-A\B~ (p + g)) = 



— e 



— e 



TUB 

rripfp 



{Ifv{Sl + Sr) + ^ Fa{Sl - Sr)^ (17) 



TUB 

in terms of the photon-helicity amplitudes Sl^r. The above relation differs from the one 
given in [23] by an overall sign, which is due to the different convention used in [23] (and 
in [51]) for the covariant derivative: = — ieQfA^ instead of = d^, + ieQfA^ 
as in this paper. In QCDF, Fax are given by QuJb/^b and induce a term (5u'^2^^ in 
a^'^^^{p~). The long-distance photon contribution to Fy^A was found to be [23] 

^SOft _ _Q Q7 ^ Q Q2 ^ Q^Q^ ^ ^SOft _ _Q Qg ^ Q Q2 = Q^Gy ■ (l8) 

with Ga + Gv = -0.24 ± 0.06 and Gy - Ga = -0.030 ± 0.015. Again, there is a relative 
sign with respect to the results in [23]. This comes from the fact that the product eF^y 
is independent of the sign convention for e, and as we have changed the overall sign of 
(17) with respect to [23], we also have to change the sign of F^y. Stated differently: 
the relative sign between F^v aiid F^v^ in [23] is wrong because of a mismatch in sign 
conventions for e in the covariant derivative. 

In order to obtain concise expressions for cl^'^r-^, it proves convenient to define one 
more hadronic quantity: 

9iR = ^ ^ (Gy ± Ga) (19) 



and correspondingly for other mesons. is 0(l/m^) as Gy + Ga has the same power 
scaling in as Ti, i.e. ~ nif^^^'^ , as one can read off from the explicit expressions in 
[49]. The difference Gy — Ga, on the other hand, is a twist-3 effect due to three-particle 
light-cone DAs of the photon and is suppressed by one more power of rrib, i.e. gR ~ 1/m^. 
This quantity will enter the CP asymmetry. Our final expressions for a^^R-^ then read: 

a^riV) - ^rr{V)\^^^^ (&^ - gl, ci^ - 0) . (20) 
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Numerically, one has (/^/fe^ = —0.3, so these corrections, despite being suppressed by 
one more power in 1/mf,, are not small numerically and larger than the known 0{as) 
corrections to QCDF from B — > 'yiu. Based on this, we feel justified in including these 
long-distance corrections in our analysis, while dropping the radiative ones of Figs. 2(a) 
and (c). 

We conclude this section by listing the numerical values of some of the annihilation 
coefficients, for central values of the input parameters, including in particular those to 
which (5i,2 contribute (with no Cabibbo suppression): 

a^r^iK*") = -0.013 - 0.001 LD , a^'l''''{K*-) = 0.004 + 0.001 LD , 
«7r"(P°) = -0-001 - 0.004 LD , a^r^iP') = 0.149 - 0.043 LD , 
a^^^''°(a;) = -0.024 + 0.003 LD . (21) 

The contribution from the long-distance photon emission is labelled "LD" (LD— 1 at the 
end). The unexpectedly small al^£^'\p^) is due to a numerical cancellation between the 
charged- current and penguin-operator contributions. Comparing these results with those 
from QCDF, Eq. (11), it is evident that WA is, as expected, largely irrelevant for the 
branching ratios, except for — > p^j. 



4 Long-Distance Contributions from Quark Loops 

In this section we calculate the soft-gluon emission from quark loops shown in Fig. 1. 
Again, these contributions are suppressed by one power of rrib with respect to a^f^'^^^, 
but they also induce a right-handed photon amplitude which is of the same order in 1/m.b 
As we shall see in the next section, this amplitude induces the time-dependent 
CP asymmetry in B ^ V^. The asymmetry is expected to be very small in the SM and 
oc mo/iTii, due to chiral suppression of the leading transition, but could be drastically 
enhanced by NP contributions. It was noticed in Refs. [17, 18] that the chiral suppression 
is relaxed by emission of a gluon from the quark loop, which is the topic of this section. 
The task is then not so much to calculate these contributions to high accuracy, but to 
exclude the possibility of large contributions to the CP asymmetry. For this reason we will 
be very generous with the theoretical uncertainties of the results obtained in this section 
— which are currently unavoidable due to the uncertainties of the relevant hadronic input 
parameters. 

Historically, soft-gluon emission from a charm loop was first considered in Rcf. [52] as 
a potentially relevant long-distance contribution to the branching ratio of -B ^ K*j, at 
about the same time as similar effects were being discussed for its inclusive counterpart 
B — > Xg'j [53, 54, 55]. It was pointed out later, in Ref. [18], that the same diagram 
also contributes dominantly to the time-dependent CP asymmetry in S° — > K*'^^. The 
size of this contribution was calculated only very recently, in Ref. [19]. The method used 
in [19] relies on the local operator-product expansion of a heavy-quark loop in inverse 
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powers of the quark mass and hence cannot be used to calculate soft-gluon emission from 
light-quark loops, which are doubly Cabibbo suppressed for 6 ^ 57 transitions, but not 
for b — > d'j. In Sec. 4.1 we will briefly review the results for heavy-quark loops and in 
Sec. 4.2 we will present a new technique for calculating light-quark loops; however, before 
we do so, we would like to fix our notation and give explicit expressions for aj^^y 

Potentially the most important contribution to the soft-gluon emission diagram in 
Fig. 1 comes from the charged-current operator with the large Wilson coefficient 
C2 ~ 1; it vanishes for QY by gauge invariance. In order to calculate the diagram, it 
proves convenient to decompose Q2 by a Fierz transformation into 

qY = ^gf + 2gf 



with = {U^U)v-a{dH)v-a. 



(22) 



The contribution of the [/-quark loop to the B — > amplitude can then be written as 



A{B V^)qv = 



Gf .(d) 
71 ^ 



C2-{V{ph{q)\2QY\B{pB)) 



C2 ■ {-ie)e*''{q) 



E 



QfJ 



dV%Vip)\T{q^^q){x)2QXmBiPB)) (23) 



a/2 



(24) 



where the minus sign comes from the sign convention for e as discussed in Sec. 2. We 
decompose Fy^ into contributions from the photon-helicity amplitudes Sl,r, Eq. (7), as 



V'^B-lu{V)P + lu{V)P 



(25) 



with 



P^e^^p<'e*Xp,q, = ]^{SL + SR), 



P = i{{e*r)*){pq) - {e*p){rj*q)} = -{S^ - Sr) . 



(26) 



In addition to Q2, the penguin operators Q 3,4,6 give a non-zero contribution to soft-gluon 
emission. Including all these contributions, and comparing (24) with (6) and (7), we 
obtain the following expression for 



C/,soft 



m,rf-^(o) 



QuC2{lu ± iu){V) + QdCs{Id ± IdW) 



q ) 



(27) 
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Here the sum over q runs over all five active quarks d, s, c, b. D denotes the down-type 
quark in the b D'-f transition. The minus sign in front of Cq is due to Furry's theorem 
according to which only the axial-vector current in the UU term in (22) contributes to 
^VB- We do not include the contribution from Q5 because its Fierz transformation changes 
the chirality of the current so that the resulting loop contribution is proportional to 
and hence helicity suppressed. In the following we distinguish between heavy {b, c) and 
light {u,d,s) quark loops. Assuming SU(3)-flavour symmetry, one has lu — Id — h, ^-nd 
ditto for lu,d,s, which causes a cancellation of these contributions in the last term in (27). 
As to be discussed below, in Sec. 4.2, we estimate the SU(3)-breaking effects to be around 
10%. 

We now turn to the calculation of l^^c, h,c and /„, 
4.1 Heavy- Quark Loops 

The calculation of lc{K*) and lc{K*) was presented in Ref. [19]; here we briefly recapitulate 
the method and present results also for If,, and for p, uj, K*, 0. 

It was first noticed in Ref. [52] that soft-gluon emission from a charm loop, Fig. 1(b), 
is suitable for an operator product expansion (OPE) in l/rric since the on-shcll photon is 
far away from the partonic threshold Atu^. The OPE reads, to leading order in l/rric [52]: 

ie" J d^xe''^-''Tc-f^c{x)2Ql ^ Qf + 0(l/m^) (28) 

where 

Qp = cpD^,{l - ^,)^gGl^D^{F"(')b 

with 

Cp = -l/(487r2m^) (29) 

and Fa/3 — i{<la^*(j — 9/3^*) for an outgoing photon. Note that here the sign of g corresponds 
to the covariant derivative D^ — d — igT"-A^ which differs from the sign convention used 
in Sec. 2, but agrees with that used as a standard in hard-perturbative QCD calculations. 
Our final results for 1^, however, are independent of the sign of g, as the matrix element 
of Qp over mesons will be expressed in terms of three-particle light-cone meson DAs 
containing an explicit factor g which refers to the same convention. 
The matrix element of can be expressed in terms of Ic, Ic as 

^VB = {V{p)\Qf\B{pb)) = h{V)P + l{V)P. (30) 

The parameters lc{K*) and lc{K*) were first calculated in Ref. [52] from three-point QCD 
sum-rules. In Ref. [19] we calculated them from LCSRs, which are more suitable for 
the problem than three-point sum rules, see the discussion in [19]. The sum rules were 
obtained for the quantities L and L, which are related to Ic and Ic by 

cpL ^l-lc, cpL = J Ic (31) 
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with cf given in Eq. (29). The LCSR for L reads [19] 



m\fB 





) 



2 



Le 



R2{u) ,(32) 



where Ri and i?2 ai'e given in terms of three-particle twist-3 and 4 DAs of the vector 
meson. Exphcit expressions are given in [19]. The sum rule for L is analogous. 

In this paper we update the values of lc{K*) and lc{K*) as determined in Ref. [19] and 
also calculate these parameters for the other vector mesons. The twist-3 and 4 parameters 
entering Ri^2 are given in Tabs. 5 and 6. The results for Ic and Ic are given in Tab. 7 ^. 
Those for lb and lb are obtained as 



Although the total uncertainties in these parameters are rather large, their contribution to 
Qj]^ is of 0(2%) at best and has an only minor impact on the branching ratios. Soft-ghion 
emission is also irrelevant for isospin asymmetries; its main impact is on the CP asymmetry 
which is small by itself, so even large uncertainties are acceptable if the aim is to rule out 
a numerically sizable time-dependent CP asymmetry in the SM. In obtaining these results 
we use the one-loop pole mass nib — (4.7 ± 0.15) GeV, fs and fs, as given in Tab. 3, the 
Borel parameter = (12 ± 3) GeV^, the continuum threshold sq — (35 ± 2) GeV^ and 
the renormalisation scale /i^ = — ± 1 GeV^ [30] . The uncertainties of the DAs are 
given in Tab. 5. Within the accuracy of the sum rules for these parameters it is impossible 
to distinguish between p° and u, so we assume them to be equal. In view of the large 
uncertainties associated with the parameters of the three-particle DAs we have adopted 
a conservative way to estimate the total uncertainty of Ic and Ic in Tab. 7 and added the 
uncertainties linearly. The uncertainties are sizable in m^, fs, Cl and M^. It is worth 
noting that the differences lc{V) — lc{V) hardly depend on the Borel parameter M^. 

Let us turn to the issue of the convergence of the 1/m.c expansion which was discussed, 
for the inclusive case, in Refs. [54, 55]. Higher-order terms in the expansion of (28) contain 
operators of type D{q ■ D)"'Gb; the expansion can be resummed with the result given in 
Ref. [55]. For inclusive decays, the relevant matrix elements are {B\bD"'Gb\B) , which 
can be estimated, on dimensional grounds, as {B\bD"'Gb\B) ^ Aqqj^ {B \bGb\ B) . The 
expansion parameter is then t = (mfeAQCD)/(4wi^) ~ 0.2, which is not power suppressed, 
but not large numerically. For t = 0.2 the effect of resummation is to enhance the leading- 
order matrix element by 15%, whereas for t — 0.4 it amounts to a 30% enhancement. We 
expect the resummation to have a similar effect in exclusive decays. We shall include the 

^The values obtained in [52] with local three-point sum rules are lc{K*) — — (1374± 250) keV and 
lc{K*) = —(1749 ±250) keV. The quoted uncertainty includes solely the variation of the Borel parameter 
and therefore probably underestimates the uncertainty. The central values arc substantially larger than 
those obtained from LCSRs, Tab. 7. It is well known that three-point sum rules are inappropriate for 
b transitions since higher-order condensate contributions grow with rUb and destroy the hierarchy of 
perturbative and non-perturbative contributions. 




(33) 
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0.65(25) 


0.35(10) 


0.26(10) 







-0.05(2) 






Table 5: Three-particle twist-3 hadronic parameters at the scale /U = 1 GeV. The parameters A 
and K are G-odd whereas the parameters and cu are G-evcn. The results for K* are updates 
of those published in [19], those for p updates of [40] and those for cp are new. We assume the 
parameters of p and lo to be equal. A full derivation of these results will be published elsewhere. 
Note that the absolute sign of all these parameters depends on the sign convention chosen for 
the strong coupling g, see text, and that K* refers to an (sq) bound state. 





p, LO, (f) 


K* 


Remarks 


Civ 


0.10 ±0.05 


0.10 ±0.05 


from [56]; no SU(3) breaking; to be updated in [57] 


Civ 


= ~Civ 


= "Civ 


ditto 







0.012(4) 


G-odd; quoted from [43] 



Table 6: Three-particle twist-4 hadronic parameters at the scale /x = 1 GeV. The same remark 
about the absolute sign and the meaning of K* applies as for twist-3 parameters. 





Ic 


Ic 


Ic Ic 


Ic + Ic 


B 


K* 


-355 ± 280 


-596 ± 520 


242 ± 370 


-952 ± 800 


B 




-382 ± 300 


-502 ± 430 


120 ± 390 


-884 ± 660 


Bs- 


^ K* 


-347 ± 260 


-342 ± 400 


-4 ± 300 


-689 ± 600 


Bs- 




-312 ± 240 


-618 ±500 


306 ± 320 


-930 ± 750 



Table 7: Soft-gluon contributions from c-quark loops in units keV. The quantities Ic and Ic arc 
defined in (24) and (25). We assume equal parameters for p and oj. li, is obtained as = Icml/ml 
and correspondingly for i^. 
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effect of truncating tlie l/nic expansion by doubling the theoretical uncertainty of our 
final result for the CP asymmetries, which depend on — lc\ the impact of l^. + If. on the 
branching ratios is small. We also would like to mention, as noted in [55] , that besides the 
derivative expansion in the gluon field there are further higher-twist contributions from 
e.g. two gluon fields. These contributions, however, are truly power suppressed and of 
order Aq^jj/m^, so we feel justified neglecting them. 

4.2 Light-Quark Loops 

For light-quark loops the photon is almost at threshold and the local OPE does not apply, 
unhke the case of heavy quarks discussed in the previous subsection. In this subsection 
we develop a method for calculating these contributions which starts from the calculation 
of Ty^{q^) from LCSRs for an off-shell photon momentum ^ 0. We then shall use 
a dispersion relation to relate the off-shell matrix element to ry^(O), which in turn can 
be expressed in terms of the wanted quantities luiV) and lu{V), Eq. (25). The starting 
point of the method developed here is similar to the one used by Khodjamirian for the 
calculation of soft-gluon contributions to i? ^ tttt [58]. In order to simplify notations, we 
drop the superscript u on the correlation function. 

A suitable correlation function for extracting ryg(O) for the weak interaction operator 
Q2 is 

(34) 

with Qi defined as in the previous subsection and pb = p + q, P = pb — k. The current 
Jb — rribbij^q is the interpolating field of the B meson with 

{B{pB)\JB\0)=mlfB. (35) 

The leading-order contribution to this correlation function, with a soft gluon, is shown in 
Fig. 3. Following Ref. [58], we have introduced an unphysical momentum k at the weak 
vertex. This additional momentum serves to avoid unphysical low-lying cuts in p^, also 
known as parasitic terms. We will choose the momentum configuration in such a way 
that k disappears when extracting ry^(O). The kinematics of the correlation function 
describes a 2-2 scattering process and therefore depends on six independent momentum 
squares. Three of those, namely 

P\{q- kf , {pI - ml) « -AjcD (36) 

are chosen to lie below their respective thresholds assuring that the correlation function 
is dominated by light-like distances and therefore suitable for a light-cone expansion. The 
other three independent variables arep^, k"^ and (f. Neglecting higher-order corrections in 
the vector-meson mass, we set = and, for simplicity, A;^ = 0. We also set q^ — {q — kY, 
which will be necessary for avoiding a subtraction constant in the dispersion relation in 
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^aaa; 



Figure 3: Leading contribution to the correlation function Ty in (34). The black square denotes 
insertion of the operator Qi with i = I, . . . ,Q. The B meson momentum is = p + q and the 
vector meson carries the momentum p. 

and also leaves only one remnant of the presence of the unphysical momentum k: 
— {Pb — kY Pb- The rationale for this choice of kinematics will become more 
transparent below. 

Inserting a complete set of hadron states, the correlation function becomes 

^v(g^p|, P') = {mlh) ^^^^ + ..., (37) 

rriB -Pb 

where the dots stand for higher states and the star on Fy^ is to remind one of the presence 
of the unphysical momentum k in P^. We can decompose the correlation function as 

Fy = 7yP + 7vP + 0(fc), 

with the projectors P and P given in Eq. (26). Additional structures in k are unphysical 
and can be dropped. Calculating the w-quark loop to twist-3 accuracy, we get 

fvmlmv f (P^ - (g - A;)^^ 



^ m-mi) ^^^^ ' ^''^ 

fymlmy r v{{q-kf-P^) 

where I = q — k + va^p and pi, = q + aip, cii = 1 — cti and therefore 

= va^P^ + (1 - vas) {q-kf , pl^aiq^^ aip\ , (40) 

where in our choice of kinematics {q — kY — > q^ in the sequel. V and A are twist-3 
three-particle DAs of the vector meson; they are discussed in detail in Ref. [40]. The 
quantities jy and 'jy also receive contributions of higher twist, which we do not include 
in this paper. The integration measure is defined as 

/ — dv daida2das6{l — a\ — a^ — aa) . (41) 

J{v,a) Jo Jo 
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Eqs. (38), (39) and (40) clearly show that the introduction of the unphysical momentum 
k avoids a low- lying cut (parasite) from the -u-quark loop in the variable p^. 

The parasitic term in q^, however, which originates from the 6-quark propagator going 
on-shell, is not absent for our choice = {q — k^. It induces a parasitic term to be added 
to (37) which is of the form 

p*P r 

{V{p)\JB\BD{q)) , _ , / d^xe'^'^-^>%BD{q)\T2d,mu^pu]{xm + 0{k) . (42) 

The matrix element on the left-hand side is just the form factor Aq{p\) ioi B ^ V 
transitions, which was calculated from LCSRs in [30] and exhibits a pole ~ l/(m^ ~ P%) 
in inducing a parasitic contribution to the first term in Eq. (37) . Before wc can proceed 
any further, wc need to determine the size of this parasitic contribution. If wc were dealing 
with a c-quark loop, we could apply a local OPE to the integral and calculate its value 
from (28) and the following estimate based on dimensional analysis: 

{Bd\D^,{1 - ^,)G^pb\0) ~ fB{PB)pt^p ■ A|cD , (43) 

where tap is an antisymmetric dimensionless tensor. The ratio of the parasitic term (42) to 
the main term in (37) is then of order AgQ^/m^ ~ 1% and neghgibly small. For the u loop, 
on the other hand, the local OPE is not apphcable and one is back to our initial problem 
of devising a method to calculate a non-local correlation function, although in this case 
a simpler one than that in Eq. (23). Given, however, the smallness of the parasitic term 
for heavy quarks ~ 1% , it is unlikely that this term is an order of magnitude larger for 
light quarks, especially in view of the numerical closeness of the light-quark and heavy- 
quark loops. For our cases of interest even a contamination at the level of 50% would not 
constitute a major problem, as these contributions to B ^ Vj are only relevant for the 
time-dependent CP asymmetry which is expected to be near zero in the SM. Our major 
aim is to confirm that this is indeed the case and to exclude large contributions from 
soft-gluon emission, but not to give a precise determination of their size. In view of this, 
even a large parasitic contamination is perfectly acceptable. 

The next step is to write (38) and (39) in terms of a dispersion relation in p^, 

Mq^pI, P') = - r ^^slv{<l\ s, P') , (44) 

in order to match them to the hadronic representation (37). The quantity l*^{q^^P'^) 
can then be obtained by applying the standard QCD sum rule techniques, namely Borel 
transformation and continuum subtraction, which yield 

m^P') = r rf«e(-^-)/^^Im,7y(?^s,n , (45) 

mBjB TT J^2 

where the star again indicates the presence of the unphysical momentum A; in = 
{Pb — kY; also note that the photon is still off-shell. Once lm.s^v{(f, P^) is known. 
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/* (g^, P^) can be analytically continued in + iO. For the 5-meson ground state 

this removes the last trace of the unphysical momentum k. The analytic continuation 
in P^ is justified because it is far above the other hadronic scales in the corresponding 
channel. This yields the physical w-quark amplitude, yet still for an off-shell photon. 

After the conceptual outline given above we will now outline how to proceed from 
the intermediate results (38) and (39). In order to obtain the imaginary parts of 7y and 
7y, it proves convenient to perform some of the integrations over the variables v and ctj 
until logarithms appear whose imaginary parts (cuts) can easily be identified. It turns 
out that the integrals over dv and dai with 0:2 = 1 — cti — as are elementary since the 
involved variables are space-like which guarantees the absence of singularities. We obtain 
an expression of the form 

7y ~ ^ dag^^^— ^{((ln[m^ - p|] - ln[m^ - a^pl - as?']) A + P2) 

(ln[-g2] - ln[-a3P' - «3g']P3 + P4)}P5 , (46) 

where Pi stands for polynomials. The poles in are integrable, i.e. removable. The 
dispersion representation in p\ is now obtained from the cuts of the logarithms. We find 



2 2 



- j J da^ ((ln[-g2] - InhaaP' - a3g'])P3 + ^4) 



P1P5 



^1 J ^-^oyv-n -11 -n '^■'^ ' ^ V(P2_g2^(^|_^2)3 

ml 

The integral over da^ is elementary and we finally obtain the imaginary part for the 
dispersion relation (44): 



-Ims7y(?^,s,P^) 

TT 



fvmlmv r 21 , 1 r 21 

(ln[-g'^J + ln[s - q^\ 



s>ml 87r2(P2 _ ^2)3(, _ ^2)5 

- ln[-m^P^ - q\s -ml- P')\ + P,)P, . (47) 

This is the expression to be used in Eq. (45). As discussed above, the momentum k 
completely disappears upon analytic continuation of P^ — > m\ + iO and we obtain the 
amphtude for an off-shell photon. The analytic continuation is rather straightforward: Z* 
acquires an imaginary part from those logarithms whose arguments depend on P^. The 
imaginary part is proportional to the mass of the u quark and originates from the quark 
going on-shell. After the analytic continuation of (47), all remnants of the unphysical 
momentum k have disappeared and we can drop the star from now on: 

W)=lu{q\ml + iQ) = ^-- r cise("^B-)/A^^ini,7y(g2,s,m^ + ^0), (48) 

m%}B TT J ml 

for q^ <^ — -^QCD- It is interesting to note that if one does not project onto the B ground 
state the analytic continuation leads to unphysical cuts in negative q^ which come from 
the fact that for higher states the unphysical momentum k is still present. 
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There remains only one step to be done, namely to put the photon on-shell, i.e. 
0. To do so, we follow the method used in Ref. [59], where the pion-photon-photon 
transition form factor F^*^t^ was estimated with one on-shell photon. Since /^(g^) is an 
analytic function in g^, it has the standard dispersion representation 



2, 1 p,Jm,/„(t) 



W) = - dt-^ (49) 

for below the cut. Potential subtraction terms spoiling the above representation are 
absent. This can be seen as follows: for very large Euclidian vahics —q^ ^ ^qcd 
one can perform a local OPE very much the same way as for c-quark loops with the 
expansion coefficient — >■ 1/q^. Using this result we have explicitly verified that 

^w(?^)) Iu{q^) ~ Another indication, although not sufficient, is that the explicit 

calculation of /* (g^, P^) does not contain a constant or polynomial terms in q^. The imag- 
inary part in q^ comes from the logarithms in (47) and the poles in l/{q^ — s). The poles 
in = m\ -\- iO are again integrable or removable. 

The perturbative or parton representation has a cut starting at 0, Eq. (49), and it is 
therefore impossible to set = because it is right below the perturbative threshold. 
The idea is then to cut out this lower part by inserting resonances that couple to the 'U7^'U 
current; = is then sufficiently below the resonances and the corresponding continuum 
threshold. We shall content ourselves with the two lowest resonances p and uj. Treating 
them as equal, we have 

,.tf)=^,iri^ (50) 

m^-q^ TT j,P t - q^ 

where 

r,P + f,P^ e;J2{^\uru\p){pV\2Q^\B) , 

pol 

and the sum runs over the polarisation of the p. It remains to determine (and fp), 
so that Eq. (50) can be used to extract /^(O). This can be achieved by applying a Borel 
transformation in the variable 5^ which yields the estimate 

p 

2r„ = - r diImZ„(i)eK-*)/A^' (51) 

and finally 

lu ^ lu{0) = - r 4^^"^'"*^^'"' ^u{t) + - r^lm Ut) . (52) 

Jo ITT'p J,P t 

The crucial point here is that for t below the continuum threshold the factor 1/t gets 
replaced by e*^™'^~*^/^^/m^. 

At this point we would also like to clarify in what respect our method to calculate 
soft-gluon emission in P — > differs from that developed in Ref. [58] for analogous con- 
tributions to the non-leptonic P — > tttt decay. In both cases the problem is a light-quark 
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~l~ 


B - 


> K* 


536 ± 70% 


635 ± 70% 


-99 ± 300 


1172 ± 70% 


B - 




827 ± 70% 


828 ± 70% 


-1 ± 300 


1655 ± 70% 


Bs- 


K* 


454 ± 70% 


572 ± 70% 


-118 ±300 


1025 ± 70% 


Bs- 


-0 


156 ± 70% 


737 ± 70% 


-581 ± 300 


893 ± 70% 



Table 8: Soft-gluon contributions from tt-quark loops in units keV. The quantities and Z„ are 
defined in (24) and (25). We assume lu{p) = lui^^) and similarly for Z^. The uncertainty for 
lu — lu is given in absolute numbers because of cancellations. In the SU(3)-flavour limit assumed 
in this calculation one has 1^ = = Ig = Iq. 

loop which is almost on-shell, the corresponding non-perturbative effects are estimated 
from light-cone sum rules and, in order to avoid parasitic terms in the correlation func- 
tion, an auxiliary momentum is introduced into the weak vertex. The distinction is that, 
in contrast to the pion, the photon is a perturbative state and therefore cannot be rep- 
resented by an interpolating current, but appears directly in the diagram with on-shell 
momenta. In order to set the photon on its mass shell we use a dispersion representation, 
Eq. (49), and estimate the truly non-perturbative part of the spectral function from the 
corresponding sum rule, Eq. (51). Moreover we have checked, by inspecting the OPE in 
the deep Euclidian, that the dispersion representation has no subtraction terms, which is 
implicitly assumed in (49). In order to assure the absence of these terms we had to set 
the two in principle independent momentum squares and (q — k)'^ equal to each other. 
This reintroduced a parasitic contribution of the form (42) which we estimated to be of 
0(1%) as compared to the main contribution. 

The sum rule (52) gives the numerical results collected in Tab. 8. We use the Borel 
parameter = (1.2 ± 0.3) GeV^ and the threshold = (1.6 ± 0.1) GeV^. Comparing 
these results with those from the c loop. Tab. 7, we see that they are roughly of the same 
size, but come with opposite sign. The smallness of (/„ — lu){p) is due to an accidental 
numerical cancellation. The uncertainties are large, which is no cause of concern, however, 
because we are only interested in the approximate size of these contributions which set 
the size of the time-dependent CP asymmetry in B ^ Vj. 

In the above, we have assumed SU(3)-fiavour symmetry which implies lu = la = h = Iq- 
We can estimate the size of SU(3)-breaking effects by taking into account that Ig couples 
to the S7^s current via the and higher resonances in (50) and requires a slightly higher 
continuum threshold Sq. This leads to a numerical difference with respect to Z„ which is 
around 5%. The effect of neglecting the quark masses is of order mq/mt, and therefore 
even smaller. We conclude that it seems unlikely that Is differs from lu,d by more than 
10%. 
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5 Phenomenological Results 



In this section we combine the different contributions to the factorisation coefficients 
"7L(i?) calculated in Sees. 2, 3 and 4 and give results for the observables in 5 — > 
transitions, namely the branching ratio, the isospin asymmetry and the time-dependent 
CP asymmetry. 



5.1 Branching Ratios 

The (non-CP-averaged) branching ratio of the b — > decay B — > is given by 



B(B - ^7) = ? Glamlml / ^ ^ 



327r4 



u 



u 



(53) 



with the isospin factors Cp±^K*,4> — 1 and CpO = -\/2- The branching ratio for the CP- 
conjugated channel B — > (b — > at parton level) is obtained by replacing A[f ^ — > 
(A[f^)*. Experimental results for B — > K*^ and i? {p,uj)'y are collected in Tab. 1. For 
Bs 07 there is only an upper bound B{Bs (j)-f) < 120 x lO'^ [28]. No experimental 
information is available for Bg K*'y. 

With the input parameters from Tab. 3 and the lifetimes given in Tab. 9 we find the 
following CP-averaged branching ratios for B — > K*'f, making explicit various sources of 
uncertainty: 

B{B- K*--f) = (53.3 ± 13.5(Ti) ± 4.8(/x) ± 1.8(^6) ± 1.9(i„,c) ± 1.3(other)) x 10"^ 



= (53.3 ± 13.5(Ti) ± 5.8) x 10"^ , 
B(S° ^ X*%) = (54.2 ± 13.2(ri) ± 6.0(/x) ± 1.8(^6) ± 1.8(Z„,c) ± 1.4(other)) x 10 



-6 



(54.2 ± 13.2(ri) ±6.7) x 10" 



(54) 



We have added all individual uncertainties in quadrature, except for that induced by the 
form factor. The uncertainty in // is that induced by the renormalisation-scale dependence, 
with = mij{mij) ± 1 GeV. The uncertainty in lu,c refers to the soft-gluon terms calculated 
in Sec. 4. "Other" sources of uncertainty include the dependence on the parameters in 
Tab. 3, on the size of LD WA contributions and the replacement of NLO by LO Wilson 
coefficients. The above results agree, within errors, with the experimental ones given in 
Tab. 1, within the large theoretical uncertainty induced by the form factor. 

As the uncertainties of all form factors in Tab. 3 are of roughly the same size, one 
might conclude that the predictions for all branching ratios will carry uncertainties similar 
to those in (54). This is, however, not the case: the accuracy of the theoretical predictions 
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Tb±/tbO 




1.530(9) ps 


1.071(9) 


0.958(39) 



Table 9: B lifetimes from HFAG [6] 



can be improved by making use of the fact that the ratio of form factors is known much 
better than the individual form factors themselves. The reason is that the values given in 
Tab. 3, which were calculated using the same method, LCSRs, and with a common set of 
input parameters, include common systematic uncertainties (dependence on fs, etc.) 
which partially cancel in the ratio. In Ref. [22] we have investigated in detail the ratio of 
the K* and p form factor and found 

^ = 1.17 ±0.09. (55) 

The uncertainty is by a factor 2 smaller than if we had calculated from the entries in 
Tab. 3; an analogous calculation for uo yields 

eu.^^p^ = 1-30 ±0.10. (56) 

The difference between and is mainly due to the difference between and see 
Tab. 3. For the Bg form factors, we also need the ratio of decay constants fsJfBd- The 
status of Jb from lattice was reviewed in Ref. [33] ; the present state-of-the-art calculations 
are unquenched with Nf = 2 + 1 active flavours [60], whose average is fss/fsi = 1-23 ± 
0.07. Again, this ratio is fully consistent with that quoted in Tab. 3, but has a smaller 
uncertainty. We then find the following ratios for Bs form factors: 



\ ' ^ 1.01 ±0.13, C^* = o r,: = 1.09 ±0.09. (57) 
^ Tf«^'^(0) rf^^^*(o) 



The uncertainty of ^x* is smaller than that of because the input parameters for K* 
and K* are the same (except for G-odd parameters like aj;) and cancel in the ratio; the 
uncertainty is dominated by that of /bJIb^- 

To benefit from this reduced theoretical uncertainty in predicting branching ratios, 
one has to calculate ratios of branching ratios, which mainly depend on and only 
mildly on Ti itself: in addition to the overall normalisation, Ti also enters hard-spectator 
interactions and power-suppressed corrections, whose size is set by hadronic quantities 
oc 1/Ti. As these corrections are subleading (in ag or l/rrib), however, a small shift 
in Ti has only very minor impact on the branching ratios. The absolute scale for the 
branching ratios is set by the CP- and isospin-averaged branching ratio with the smallest 
experimental uncertainty, i.e. B — > K*^; from Tab. 1, one finds: 



B{B K*-f) = - { B{B^ K*^-t) + — fi(B" ^ i^*%) ) = (41.6 ± 1.7) x 10"" . (58) 

2 1 Tgo 
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That is, we obtain a theoretical prediction for B{B ^ V-y) as 

B{B V-f) 



B{B^Vj)\ 



th 



B{B K*-f) 



B{B K*-f) 



th 



exp 



(59) 



where 



Jth 



depends mainly on and only in subleading terms on the individual form 
factors Tf'*^* and Tf~'^. It is obvious that, except for these subleading terms, this 
procedure is equivalent to extracting an effective form factor Tf^^*{0)\^^ from B — > 
K*'-f and using Tf^'^{0)\^^ = Tf^^*{0)\^g/^v for calculating the branching ratios for 
B V^. From (58) we find 



T^-^^' (0) l^g = 0.267 ± 0.017(th) ± 0.006(exp) = 0.267 ± 0.018 , 



(60) 



where the theoretical uncertainty follows from the second uncertainty given in (54). 
Eqs. (55), (56) and (57) then yield 



rf^^(o) 



*K* 



(0) 



eff 



cff 



0.228 ±0.023, 
0.245 ± 0.024 , 



"""(O) leff = 0.205 ± 0.021 , 



7f^^'^(0) 



= 0.260 ± 0.036 . 



eff 



(61) 



Note that all effective form factors agree, within errors, with the results from LCSRs 
given in Tab. 3, which confirms the results obtained from this method; the crucial point, 
however, is that the uncertainties are reduced by a factor of 2 (except for Tf"'^'^). We 
would like to stress that the motivation for this procedure is to achieve a reduction of 
the theoretical uncertainty of the predicted branching fractions in 5 — > (p, cu)^ and Bg 
decays. The effective form factors do not constitute a new and independent theoretical 
determination, but are derived from the experimental results for B K*y under the 
following assumptions: 

• there is no NP in S — > K*^-^ 

• QCDF is vahd with no systematic uncertainties; 

• LCSRs can reliably predict the ratio of form factors at zero momentum transfer. 
Prom (53) and (61), we then predict the following CP-averaged branching ratios: 



B{B- - 


^ P l) = 


(1.16 ±0.22(Ti) ±0.13) X 10-^ 


F(S°- 




(0.55 ± O.ll(ri) ± 0.07) X 10-^ , 


F(S° 




(0.44 ± 0.09(ri) ± 0.05) X 10"^ 


B{Bs- 


. X*7) = 


(1.26 ± 0.25(Ti) ± 0.18) x 10"^ , 


B{Bs 


<t>l) = 


(39.4 ± 10.7(ri) ± 5.3) X 10"^ , 



(62) 



^Which is motivated by the results from inclusive B Xg'f decays [2]. 
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QCDF 


+ WA (no LD) 


+ WA (incl. LD) 


+ soft gluons 


B- - 


P 1 


1.05 


1.17 


1.11 


1.16 






0.49 


0.53 


0.53 


0.55 


5°- 


a;7 


0.40 


0.42 


0.42 


0.44 


B- - 


ir*-7 


39.7 


38.4 


38.3 


39.4 


5°- 




37.1 


39.7 


39.9 


41.0 


Bl- 


i?*07 


1.12 


1.22 


1.23 


1.26 




07 


34.6 


38.2 


38.3 


39.4 



Table 10: Individual contributions to CP-averaged branching ratios, using effective form factors 
and central values of all other input parameters given in Tab. 3 (in particular 7 = 53°). LD 
stands for the long-distance photon-emission contribution to WA. Each column labelled "-|-X" 
includes the contributions listed in the previous column plus the contribution induced by X. The 
entries in the last column are our total central values. 



where the first uncertainty is induced by the effective form factors and the second includes 
the variation of all inputs from Tab. 3 except for the angle 7 of the UT, which is fixed at 
7 = 53°. The total uncertainty in each channel is ~ 20%, except for B^ ^7, where it is 
30%. The results for p and lo agree very well with those of BaBar, Tab. 1, but less so with 
the Belle results, although present experimental and theoretical uncertainties preclude 
a firm conclusion. Our prediction for Bg — > ^7 is well below the current experimental 
bound 120 x 10~^ [28]. A branching ratio of the size given in (62) implies that O(IO^) 
Bg 07 events will be seen within the first few years of the LHC. In Tab. 10 we detail 
the contributions of individual terms to the branching ratios. In all cases B is dominated 
by the QCDF contribution, with WA most relevant for B^ P I- This is expected as 
WA enters with the large Wilson coefficient C2 ~ 1. The effect is extenuated by long- 
distance (LD) photon emission, which itself is compensated by soft-gluon emission. The 
other channels follow a similar pattern, although the size of the effects is smaller. 

Let us now turn to the determination of CKM parameters from the branching ratios. 
In this context, two particularly interesting observables are 

B{B K*'y) ' ' B{B^ K*'y) ' ^ ^ 

given in terms of the CP- and isospin-averaged branching ratios of B ^ (p, uj)^f and 
B p7, respectively, and B K*'~f decays, see Tab. 1. Rp/uj has been measured by both 
BaBar and Belle [4, 5], a first value of Rp has been given by BaBar [4]. The experimental 
determinations actually assume exact isospin symmetry, i.e. r{B^ — > p^7) = 2r{B^ 
p°7), and also r(S° p°7) = T{B^ — > lo^); as we shall discuss in the next subsection, 
these relations are not exact, and the symmetry-breaking corrections can be calculated. 
Hence, the present experimental results for Rp/^ are theory-contaminated. As the isospin 
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asymmetry between the charged and neutral p decay rates turns out to be smaller than 
the asymmetry between and c<j, it would actually be preferable, from an experimental 
point of view, to drop the uj channel and measure Rp instead of Rp/^ji as done in the 
most recent BaBar analysis on that topic [4]. We will give numerical results and theory 
uncertainties for both Rp/^j and Rp. 

One parametrisation of Rp/^j often quoted, in particular in experimental papers, is 



R 



Vtd 

Vts 



1 



1 — m\*/m\ 



So 



[1 + Ai?], 



(64) 



with Ai? = 0.1 ±0.1 [9] and again assuming isospin symmetry for p and uj. This parametri- 
sation creates the impression that Ai? is a quantity completely unrelated to and with a 
fixed value independent of |\4d/Vts|. We would like to point out here that this impression 
is wrong: AR contains both QCD (factorisable and non-factorisable) effects and such 
from weak interactions. In Ref. [22], we have expressed AR in terms of the factorisation 
coefficients a^^, assuming isospin symmetry for p° and lu, as 



1 + AR 



1 + Re {6a± + Sgq) 



Rb cos 7 



1 - 2RbCos-f + Rl 



2Rb cos 7 + 



(65) 



with 5ao,± = '27l(p'''^)/'^7l(p°'^) ~ 1- Here 7 is one of the angles of the UT (7 = argl^*^ 
in the standard Wolfenstein parametrisation of the CKM matrix) and R^ one of its sides: 



Rh 













Vcb 



Eq. (65) shows exphcitly that AR depends both on QCD (5a±,o) and CKM parameters 
(i?6,7). The point we would like to make is that the calculation of AR requires input 
values for Rf, and 7. Once these parameters (and the Wolfenstein parameter A) are fixed, 
however, iVtd/Vtsj is also fixed and given by 



Vtd 

Vts 



X^l-2RbCos-f + Rl 



1 + 2(1 



(66) 



Hence, as and {Rb, 7) are not independent of each other, it is impossibleto extract 

from (64) with a fixed value of AR. We hasten to add that our arguments rely 
on the unitarity of the CKM matrix, and its well-known consequence, the existence of the 
UT. The unitarity of the CKM matrix is, however, already hard-wired into the effective 
Hamiltonian (2); without it, the theory would look quite different because of the absence 
of the GIM mechanism, as mentioned in Sec. 2, Eq. (3). Stated differently: as long as 
Eq. (2) is adopted as the relevant effective Hamiltonian for b — > D7 transitions, unitarity 
of the CKM matrix is implied. Obviously, the unitarity of the CKM matrix is subject 
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to experimental scrutiny, but any test of it has to involve the comparison of different 
measurements described within the same framework (by the same effective Hamiltonian), 
while a mixture of different frameworks (unitary vs. non-unitary CKM matrix) within 
one observable, like Rp/uj, does not make any sense. 

Of course Rp/uj and Rp of (63) can be used in a meaningful way to extract information 
about CKM parameters, but in order to do so one has to settle for a set of truly indepen- 
dent parameters. Based on (66), one can exchange, say, 7 for |l^d/14s| So we can either 
consider Ry as a function of the CKM parameters Rh and 7 (let us call this the 7 set of 
parameters) or as a function of Rh and (to be called the \Vtx\ set). Using the 7 

set, a measurement of Ry{'~f,Ri,) allows a determination of 7, whereas Rv{\Vtd/Vts\, Rb) 
allows the determination of iV^d/Vtsl. In either case, the simple quadratic relation (64) 
between Ry and |V((i/\4s| becomes more complicated. In Figs. 4 and 5 we plot the re- 
sulting values of |Vj<i/VtsP and 7, respectively, as a function of Ry. Although the curve 
in Fig. 4(a) looks like a straight hne, as naively expected from (64), this is not exactly 
the case, because of the dependence of AR on jl^d/T^sl- In Fig- 4(b) we plot AR for the 
\Vtx\ set of parameters. The dependence of AR on is rather strong. Apparently 

indeed AR = 0.1 ± 0.1 in the expected range 0.16 < |V4(i/Vis| < 0.24, but this estimate 
does not reflect the true theoretical uncertainty which is indicated by the dashed lines in 
the figure. 

It is now basically a matter of choice whether to use Rp/uj to determine |Vt(i/Vts| or 7. 
Once one of these parameters is known, the other one follows from Eq. (66). In Fig. 5 we 
plot 7 as a function of Rp/a;, together with the theoretical uncertainties. In Fig. 6 we also 
compare the central values of Rp/^, with those of Rp, as a function of |T4d/T4s|- Although 
the difference is small, Rp is expected to be larger than Rp/^. In order to facilitate the 
extraction of |V^(i/V^s| (or 7) from measurements of Rp/ui or Rp, Tabs. 11 and 12 contain 
explicit values for the theoretical uncertainties for representative values of Rp/^j and Rp. 
The uncertainty induced by ^p is dominant. As discussed in Rcf. [22], a reduction of this 
uncertainty would require a reduction of the uncertainty of the transverse decay constants 
fy of p and K*. With the most recent results from BaBar, Rp/^ = 0.030 ± 0.006 [4], and 



BaBar: 
Belle: 



Rp/uj 






Vts 




Vtd 




Vts 





7 = (61.0l}|^(exp)l«|(th))°, 

7 = (65.7l^S:?(exp)l|^(th))°. (67) 



These numbers compare well with the Belle result [32] from tree-level processes, 7 = 
(53 ± 20)°, quoted in Tab. 3, and results from global fits [37]. We also would like to 
point out that the above determination of 7 is actually a determination of cos 7, via 
Eq. (66), and implies, in principle, a twofold degeneracy 7 <-> 27r — 7. This is in contrast 
to the determination from B — > D^*^K^*^ in [32], which carries a twofold degeneracy 



^Strictly speaking, (66) only fixes cos 7 as function of |V(d/Vts|, leaving a twofold degeneracy of 7. 
Eq. (65), however, only depends on C0S7, so that indeed one can unambiguously replace 7 by |Vtci/^s|- 
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Figure 4: Left panel: iV^d/Vj^p as function of Rp/^J, Eq. (63), in the \Vtx\ basis, see text. Solid 
line: central values. Dash-dotted lines: theoretical uncertainty induced by = 1.17 it 0.09, 
(55). Dashed lines: other theoretical uncertainties, including those induced by \Vub\^ \Vcb\ arid 
the hadronic parameters of Tab. 3. Right panel: A/2 from Eq. (65) as function of | Vtd/Vts| for the 
\Vtx\ set of CKM parameters. Solid line: central values. Dashed lines: theoretical uncertainty. 




Figure 5: The UTangle 7 as function of R^j^ in the 7 set of CKM parameters. Solid lines: 
central values of input parameters. Dash-dotted lines: theoretical uncertainty induced by = 
1.17 ± 0.09. Dashed lines: other theoretical uncertainties. 
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\Vtd/Vt,\ 




Mother th 


7 


\. 


Mother th 


0.026 


0.183 


±0.012 


±0.007 


50.8 


+7.5 
-8.2 


±5.8 


0.028 


0.191 


+0.012 
-0.013 


±0.006 


56.0 


+7.7 
-8.3 


±4.7 


0.030 


0.199 


±0.013 


±0.006 


61.0 


+7.9 
-8.4 


±4.0 


0.032 


0.207 


+0.013 
-0.014 


±0.006 


65.7 


+8.1 
-8.5 


±3.6 


0.034 


0.214 


±0.014 


±0.006 


70.2 


+8.4 
-8.8 


±3.5 


0.036 


0.221 


+0.014 
-0.015 


±0.006 


74.5 


+8.8 
-9.0 


±3.7 



Table 11: Central values and uncertainties of |Vt,//Vta| and 7 extracted from representative 
values of Rp/^, Eq. (63). A^^ is the uncertainty induced by ^p, Eq. (55), and Aotherth that by 
other input parameters, including and \Vuh\. 



Rp 


\Vu/Vt,\ 




Aothcr th 


7 


\. 


Aothcr th 


0.028 


0.186 


±0.016 


±0.005 


52.4 


+9.9 
-10.3 


±5.0 


0.030 


0.193 


±0.016 


±0.005 


57.4 


+10.2 

-10.3 


±3.9 


0.032 


0.201 


±0.017 


±0.005 


62.0 


±10.5 


±3.1 


0.034 


0.208 


±0.017 


±0.004 


66.4 


+10.8 
-10.7 


±2.7 


0.036 


0.215 


±0.018 


±0.004 


70.7 


+11.3 
-11.0 


±2.5 



Table 12: Ditto for Rp. A^^ is larger than in Tab. 11 because of the increased weight oi B ^ pj 
in the isospin average; Aotherth is smaller because does not enter. 



7 ^ TT + 7. Obviously these two determinations taken together remove the degeneracy 
and select 7 55° < 180°. If 7 fa 55° + 180° instead, one would have \Vtd/Vts\ ^ 0.29 
from (66), which is definitely ruled out by data. Hence, the result (67) confirms the SM 
interpretation of 7 from the tree-level CP asymmetries in B ^ D^*^K^*\ 

We would like to close this subsection by making explicit the dependence of the three 
B — > (p, a;)7 branching ratios on 7. In Fig. 7 we plot these branching ratios, for central 
values of the input parameters, as functions of 7. We also indicate the present experi- 
mental results from BaBar [4], Tab. 1, within their la uncertainty. 

5.2 Isospin Asymmetries 

The asymmetries are given by 
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Figure 7: CP-averaged branching ratios of B ^ (p, uj)'y as function of 7, using the effective 
form factors and central values of other input parameters, (a): — > p^7, (b): B^ — > ^^^7, 
(c): B^ ujj. The boxes indicate the la experimental results from BaBar [4], Tab. 1. Note 
that the resulting value of 7 from the average of all three channels is 7 = (61.0^Jg'Q(exp)^g'2)°, 
Eq. (67). 



2r(^° ^ p°7) 
r(5± ^ p±7) 

AAK*) = =^ J ^ ; (70) 

the partial decay rates are CP-averaged; Aj{p), Ai{K*) are isospin asymmetries. 

Let us first discuss A{p,uj) and Aj{p) which are relevant for the experimental deter- 
mination of B{B — >■ (p, cc;)7), which in turn is used for the determination of |Vtrf/Vts| (or 
7), see Sec. 5.1. The present experimental statistics for b ^ dj transitions is rather low, 
so the experimental value of B{B {p,Lj)'y) is obtained under the explicit assumption 
of perfect symmetry, i.e. r{B^ P^l) = 2r(i?° p'^7) = 2T{B'^ ^l)- In reality, 
the symmetry between p° and uj is broken by different values of the form factors, and 
isospin symmetry between neutral and charged p is broken by photon emission from the 
spectator quark, the dominant mechanism of which is WA, as discussed in Sec. 3. From 
the formulae for individual branching ratios, Eq. (53), and the various contributions to 
the factorisation coefficients clj^ji^ collected in Sees. 2, 3 and 4, we find 

A(p,cj) = -0.20±0.09(th). (71) 
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40° 


50° 


60° 


70° 


Mp) 


-(5.3 ±6.9)% 


(0.4 ±5.3)% 


(5.7 ±3.9)% 


(10.5 ±2.7)% 



Table 13: Isospin asymmetry Ai{p), Eq. (69), for different values of 7. 

The uncertainty is dominated by that of the form factor ratio T^^^ [Q) / T^~^ ^ {Q) = 0.90 ± 
0.05.^ The dependence on all other input parameters is marginal. Aj{p), on the other 
hand, is very sensitive to 7, whereas the form factors drop out. It is driven by the WA 
contribution and, in the QCDF framework, vanishes if WA is set to zero. In Fig. 8(a) 
we plot Ai{p) as function of 7, including the theoretical uncertainties. As suggested 
by the findings of Ref . [48] , these results are not expected to change considerably upon 
inclusion of the non-factorisable radiative corrections of Fig. 2(c). In Tab. 13, we give the 
corresponding results for several values of 7, together with the theoretical uncertainty. 
Our result agrees very well with that obtained by the BaBar collaboration: ^/(p)BaBar = 
0.56 ±0.66 [4]. 

Aj{K*) was first discussed in Ref. [20], including power-suppressed 0{as) corrections 
which unfortunately violate QCDF, i.e. are divergent. It is for this reason that we decide 
to drop these corrections and include only leading-order terms in a^. We then find 

Ai{K*) = (5.4 ± 1.0(/x) ± 0.6(NLO ^ LO) ± 0.6(/b) ± 0.6(other))% 

= (5.4 ±1.4)%, (72) 

where NLO LO denotes the uncertainty induced by switching from NLO to LO ac- 
curacy in the Wilson coefficients and "other" summarises all other sources of theoret- 
ical uncertainty. As can be inferred from the entries in Tab. 1, the present experi- 
mental result is Aj{K*)exp = (3.2 ± 4.1)%. In Ref. [20], Kagan and Neubert pointed 
out that Ai{K*) is very sensitive to the values of the Wilson coefficients C^g^ in the 
combination ae = Cf^^ + In the SM, varying the renormahsation scale as 

fjL — mb{mb) ± 1 GeV and switching between LO and NLO accuracy for the Wilson coeffi- 
cients, one has ae — — 0.039 ± 0.008, which actually induces the bulk of the uncertainty in 
(72). In Fig. 8(b) we plot Ai{K*) as function of ae/af^, with af^ = -0.039. The figure 
clearly indicates that, although there is presently no discrepancy between theoretical pre- 
diction and experimental result, a reduction of the experimental uncertainty of Ai{K*) 
may well reveal some footprints of NP in this observable. 

^Note that this rcsuh is dominated by the ratio of decay constants given in Tab. 3 and discussed in 
the appendix. The experimental results entering these averages have a large spread which may cast a 
shadow of doubt on the averaged final branching ratios for — > e+e~ quoted by PDG [28]. 
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Figure 8: Left panel: isospin asymmetry Aj{p), Eq. (69), as function of the UTangle 7. Solid 
line: central values of input parameters; dashed lines: theoretical uncertainty. Right panel: 
Aj{K*), Eq. (70), in percent, as function of the ratio r = ae/og^ of the combination of penguin 
Wilson coefficients = Cq + C5/3. Solid line: central value of input parameters, dashed lines: 
theoretical uncertainty. The box indicates the present experimental uncertainty and the straight 
black lines the theory uncertainty in r. 

5.3 CP Asymmetries 

The time-dependent CP asymmetry in V^'j is given by 

r(50(t) ^ V7) - r(50(t) ^ v-f) 



Acp{t) 



^(V"7) sin(AmBt) - C{V-f) cos(AmBt) 



(73) 



where we have neglected the width difference AT of the two neutral B mesons. This 
approximation is well justified for Bd, but less so for Bs. Although the above formula 
can easily be adapted to non-zero AF^, we refrain from doing so: the whole point in 
calculating the CP asymmetry is not so much to give precise predictions for 5* and C, but 
rather to exclude the possibility of large corrections to the naive expectation S ~ niD/mb. 
With this is mind, small corrections from a non-zero AF^ are irrelevant. 

Let us briefly recall the reason for the expected smallness of S. In the process b —>■ D'j, 
in the SM, the emitted photon is predominantly left-handed in b, and right-handed in b 
decays. This is due to the fact that the dominant contribution to the amplitude comes 
from the chiral-odd dipole operator Qj, Eq. (4). As only left-handed quarks participate in 
the weak interaction, an effective operator of this type necessitates, in the SM, a helicity 
flip on one of the external quark lines, which results in a factor rrif, (and a left-handed 
photon) in bR Dl'Jl and a factor niD (and a right-handed photon) in 6^ Drj^. 
Hence, the emission of right-handed photons is suppressed by a factor rriD/rnb, which 
leads to the QCDF prediction (13) for Oy^. 

The interesting point is not the smallness of the CP asymmetry per se, but the fact that 
the helicity suppression can easily be alleviated in a large number of NP scenarios where 
the spin flip occurs on an internal line, resulting in a factor irii/mh instead of rriD/mh. A 
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prime example is left-right symmetric models [61], whose impact on the photon polari- 
sation was discussed in Refs. [16, 17, 18]. These models also come in a supersymmetric 
version whose effect on 6 — >■ 57 was investigated in Ref. [62]. Supersymmetry with no 
left-right symmetry can also provide large contributions to 6 — > D^r, see Ref. [63] for 
recent studies. Other potential sources of large effects are warped extra dimensions [64] 
or anomalous right-handed top couplings [65] . Unless the amplitude for h — > D^r is of the 
same order as the SM prediction for b — > -D7L, or the enhancement of ft — > D^r goes along 
with a suppression oi h ^ D^l, the impact on the branching ratio is small, as the two 
helicity amplitudes add incoherently. This implies there can be a substantial contribution 
of NP to b ^ D'-f escaping detection when only branching ratios are measured. 

Although the photon helicity is, in principle, an observable, it is very difficult to 
measure directly. It can, however, be accessed indirectly, in the time-dependent CP 
asymmetry in S° — > V^, which relies on the interference of both left- and right-helicity 
amplitudes and vanishes if one of them is absent. In terms of the left- and right-handed 
photon amplitudes of Eq. (5) one has 
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(74) 

Here q/p is given in terms of the B^-B^ mixing matrix M12, in the standard convention 
for the parametrisation of the CKM matrix, by 



q M- 



'12 



P V M12 

with, in the Wolfenstein parametrisation of the CKM matrix. 



cb 



sin 7. (75) 



This method of accessing the right-handed photon amplitude via S{V'-f) was first 
suggested in Ref. [16] and later discussed in more detail in Refs. [17, 18]. The direct CP 
asymmetry C{Vj) is less sensitive to Ar, but very sensitive to the strong phase of Al and 
vanishes if the radiative corrections to a^f^'^^^, Eq. (11), are neglected. As the accuracy of 
the prediction of strong phases in QCDF is subject to discussion, and in any case C{V'-f) 
is less sensitive to NP than S{V^), we shall not consider direct CP asymmetries in this 
paper. S{V^) is rather special in the sense that usually NP modifies the SM predictions 
for time-dependent CP asymmetries by affecting the mixing phase (as in Bg J/'il^<P, 
see for instance Ref. [66]), introducing new weak phases or moderately changing the size 
of the decay amplitudes which, in the absence of precise calculational tools, makes it 
difficult to trace its impact. In contrast, the time-dependent CP asymmetry in B^ V'y 
is very small in the SM, irrespective of hadronic uncertainties, and NP manifests itself by 
relieving this suppression. The smallness of the asymmetry in the SM, and the possibility 
of large effects from NP, makes it one of the prime candidates for a so-called "null test" 
of the SM, as recently advertised in Ref. [67]. 
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The fly in the ointment, however, is that in addition to the hehcity-suppressed contri- 
bution from Qi, Ar also receives contributions from the parton process b D'-fg, which 
come without a helicity-suppression factor [17, 18]. These contributions are dominated by 
soft-gluon and long-distance photon emission in weak annihilation and are also included 
in a^^, Eq. (9). In Ref. [18] it was inferred from a dimensional estimate that these contri- 
butions could be as large as ~ 10%, but a recent explicit calculation of the contribution 
of Q2 to S{K*'y) has shown that their true size is much smaller [19]. In this paper, we 
extend the calculation of [19] to all V^'j channels and include the effects from all 

four-quark operators in the effective Hamiltonian (2) and also the contribution from weak 
annihilation. 

With Al,r and Al.r as given in (6) wc can calculate 5* directly from (74) and obtain, 
making explicit the contributions from different sources:^ 



5(P7) = 


(^O^ + ^O^+^O^ ± 1.6)% = 


(0.2 + 1.6)%, 




mo/frib LD WA soft g 






(0.01 - 0.08 + 0.22 ±1.7)% 


(0.1 + 1.7)%, 




-(2.9-0 + 0.6±1.6)% 


-(2.3 + 1.6)%, 


S{K*^) = 


(0.12 + 0.03 + 0.11 ± 1.3)% 


(0.3 + 1.3)%, 


5(07) = 


(0 + + 5.3 + 8.2) X 10-2% 


(0.1 + 0.1)%. 



Including only the helicity-suppressed contribution, one expects, for B — > K*^^ neglecting 
the doubly Cabibbo suppressed amplitude in Xu \ see Eq. (2), 

S{K*^) gi,,„, = -2 sni cpa ^ -2.7% . (77) 

For — > 07, one expects the CP asymmetry to vanish if the decay amplitude is propor- 
tional to Xf \ which, at tree level, precludes any contributions of type sm{(f) 3)171 s/mb and 
also any contribution from WA.^ 

The actual results in (76) disagree with the above expectations because of the contri- 
butions from soft-gluon emission, which enter a^^, and, for S{(f)^), because the soft-gluon 
emission from quark loops is different for u and c loops, see Sec. 4, so that a^p^ 7^ a"^ and 
hence Ar (Al) is not proportional to A^^'' ((Aj*"*)*). Note that a substantial enhancement 
of 5(07) by NP requires not only an enhancement of \Ar\ (and |v4z,|), but also the pres- 
ence of a large phase in (74); this could be either a large Bg mixing phase which will also 
manifest itself in a sizable CP violation in, for instance, Bg — > J/'4>(t>, see Ref. [66]; or it 
could be a new weak phase in Ar (and Al)] or it could be a non-zero strong phase in 

^These results are obtained using LO Wilson coefficients. The difference between LO and NLO results 

is marginal. 

^This is because the mixing angle 4>s is given by arg[(A|''')^], Eq. (75), and the interference of ampli- 
tudes in (74) also yields a factor (Aj"-*)^, if the individual amplitudes are proportional to x\^^ or (Aj*-*)*, 
respectively; this is indeed the case for the helicity-suppressed term ms/mi, induced by the operator Q7, 
Eq. (4), and the WA contributions to a^jiiip), Eqs. (15), (20), so that the phases cancel in (74). 
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one of the coefficients. Based on the calculation in Sec. 4.2 we do not see much scope 
for a large phase in a"^ (whose contribution is, in addition, doubly Cabibbo suppressed), 
but the situation could be different for a^'^"^*, where we only included the leading-order 
term in a l/rric expansion, which does not carry a complex phase, see Sec. 4.1. It is 
not excluded that a resummation of higher-order terms in this expansion will generate a 
non- negligible strong phase — which is not really relevant for our results in Eq. (76), but 
could be relevant for the interpretation of any NP to be found in that observable. For 
S'(i^*7), on the other hand, no new phases are required, and any enhancement of \Ar\ 
(and \Al\) by NP will result in a larger value of S{K*'~f). 

For all S except S{K*'j), the uncertainty is entirely dominated by that of the soft- 
gluon emission terms lu,c — ^u,c, whose uncertainties we have doubled with respect to 
those given in Sec. 4. The smallness of S{{p,uj)'~f) is due to the fact that the hclicity 
factor is given by md/mb (we use mu,d/ms = 1/24.4 from chiral perturbation theory). 
For K*, the suppression from the small mixing angle is reheved by the fact that both 
weak amplitudes in X^^ contribute, with different strength, so that the CP asymmetry is 
comparable with that of p and u. Despite the generous uncertainties, it is obvious that 
none of these CP symmetries is larger than 4% in the SM, which makes these observables 
very interesting for NP searches. The present experimental result from the B factories, 
S{K*^) = —0.28 ± 0.26 [6], certainly encourages the hope that NP may manifest itself 
in that observable. While a measurement of the h ^ d CP asymmetries is probably very 
difficult even at a super-flavour factory, S{K*'y) is a promising observable for B factories 
[24], but not for the LHC.^ Bs 0(^ X+X~)7, on the other hand, will be studied in 
detail at the LHC, and in particular at LHCb, and any largely enhanced value of 5'(07) 
will be measured within the flrst years of running. 



6 Summary and Conclusions 

In this paper we have presented a comprehensive study of the observables in B ^ 
decays, namely branching ratios, isospin and CP asymmetries, for all Bs and B^^d tran- 
sitions,^ including the most recent results on form factors from QCD sum rules on the 
light-cone and hadronic parameters describing twist-2 and 3 two- and three-particle light- 
cone distribution amplitudes of vector mesons. Our study is based on QCD factorisation 
[9, 10, 11, 12, 13, 14], but goes beyond it by including power-suppressed non-factorisable 
corrections from long-distance photon emission and soft-gluon emission from quark loops 
which are also calculated from hght-cone sum rules. In Sec. 4.2 we have devised a method 
for calculating such soft-gluon emission from a light-quark loop for an on-shcll photon, 
building on the calculation of related effects in ^ tttt [58]. The main idea is to csl- 

^K* has to be traced via its decay into a CP eigenstate, i.e. Kstt^. Neutrals in the final state are not 
really LHC's favourites. 

^We have not included pure annihilation decays, for instance —> 07, as their SM branching ratios 
are tiny, O(10~^^), and sensitive to higher-order effects in the electromagnetic interaction, which are not 
considered in this paper, but for instance in Ref. [68]. 
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culate the loop for an off-shell photon and then use a dispersion representation to relate 
it to the on-shell amplitude. For phenomenology, light-quark loops are only relevant for 
b ^ d transitions, as otherwise they are Cabibbo suppressed or come with small Wilson 
coefficients. Our estimates may be of interest also for inclusive 6 — > (^7 transitions, where 
an interplay between exclusive and inclusive effects could take place similar to that for 
b^s-f [53, 54, 55]. 

Our main results are given in Sec. 5. We find that the theoretical uncertainty of the 
branching ratios gets reduced by exploiting the fact that ratios of form factors from QCD 
sum rules on the light-cone are known with better accuracy than the form factors them- 
selves. This allows us to predict the branching ratios of all B — ^ transitions with 
~ 20% theoretical uncertainty (except for fi, ^7 which comes with a ~ 30% uncer- 
tainty), based on the experimental input from B — >■ K*^. The effect of power corrections 
beyond QCD factorisation is non-negligible for all decay channels, although in some chan- 
nels the net corrections nearly cancel. We have determined iVt^/VtsI and, equivalently, 7, 
from the most recent BaBar [4] and Belle [5] results for B{B {p,u;)^)/B{B — > K*^) as 



BaBar: 
Belle: 



td 



Vts 

Vtd 



ts 



0.199lHi(exp)±0.014(th) ^ 7 = (61.0li|^(exp)isi(th))° , 
0.207tHi(exp)iHl5(th) - 7 = (65.7t^o':?(exp)l^i(th))° . 



As the relation (66) between |Vi(i/Vts| and 7 relies on cos 7, these results have a twofold 
degeneracy 7 ^ —7. Taken together with the tree-level CP asymmetries in i? — >■ D^*^K^*\ 
for instance 7 = (53 ± 20)° from Belle [32], which comes with the discrete ambiguity 
7 <-> 7 -|- TT, our result removes the ambiguity and confirms that 7 < 180° as predicted in 
the SM. 

As for the isospin asymmetries, we find a non-zero asymmetry for the and u chan- 
nel which is driven by the difference of the corresponding form factors. The asymmetry 
between the neutral and the charged p channel, on the other hand, is very sensitive to 7, 
neglected radiative corrections and hadronic input parameters, which precludes a precise 
statement about its size. The isospin asymmetry in i? — > K*'y depends only mildly on the 
input parameters, but is sensitive to the contribution of the penguin operators Qs.e- The 
sign of the asymmetry is predicted unambiguously. Although the present experimental 
uncertainty of the asymmetry is too large to allow any definite conclusion, any reduction 
could be translated into a constraint on NP contributions to the Wilson coefficients of 
these operators. 

The time-dependent CP asymmetry S{V^) in B'^ V^'-f is sensitive to the photon 
polarisation amplitudes and is power-suppressed in the SM. The contribution of Qj is 
helicity suppressed; the contributions of other operators enter via the parton process 
b — > Djg with no helicity suppression, but are also found to be small. The largest CP 
asymmetry Ri —2% is expected for B — > K*^, whereas all other CP asymmetries are below 
the 1% level. Any value significantly different from zero, measured either at the LHC or 
a future flavour factory, will constitute an unequivocal signal for NP with non-standard 
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flavour- changing interactions. 

We also would like to discuss other results for B V'~f available in the literature. 
Obviously, there are earlier results from SCET [8] and QCD factorisation, Refs. [11, 12, 13, 
14], with which we agree apart from the effects of the new non-factorisable contributions 
calculated in this paper and/or updated hadronic input. A variant of QCD factorisation 
has been advocated and pursued by Ali and Parkhomenko (AP) [9, 10]. Another approach 
is that of perturbative QCD factorisation (pQCD), which has been apphed to S — > 
in Ref. [15]. Most observables discussed in this paper, branching ratios, isospin and CP 
asymmetries, have been calculated in both approaches, for B {K* , p,uj)'~f, and we 
shall compare the corresponding results to ours in turn. As for the branching ratios, 
it is evident from Eq. (53) that the predictions depend primarily on the form factor Ti 
and only to a lesser extent on the specific implementation of QCD factorisation. For 
this reason, as AP use the same form factors as we, namely our predictions from QCD 
sum rules on the hght-cone [30], their results in their latest update Ref. [10] are very 
close to ours. The branching ratios obtained in pQCD, on the other hand, are by more 
than a factor of two larger than ours. This discrepancy is very likely to be caused by 
larger values of their form factors, calculated within the same formalism; a more detailed 
comparison is, however, difficult because Ref. [15] docs not give any explicit numbers for 
the Ti. Turning to isospin asymmetries, AP obtain approximately the same asymmetry 
A{p,u!), Eq. (71), between and uj as we do, for the same reason as above. The isospin 
asymmetry between the neutral and the charged p, Eq. (69), is more delicate, driven by 
weak annihilation contributions and very sensitive to the angle 7, see Tab. 13. Our value, 
(5.7 ± 3.9)% for 7 = 60°, disagrees with that given by AP, Ai{p) = -(2.8 ± 2.0)% for 
the same angle. A hkely reason is the smaller size of the weak annihilation amplitude 
obtained in Ref. [49], on which AP rely for that contribution, as compared to the QCD 
factorisation result. Indeed, reducing the size of the weak annihilation contribution in the 
B V'-f amplitude, our results move closer to those of AP. Ref. [15], on the other hand, 
obtains Aj{p) = (5.7 dz 6.0)% for 7 ^ 60°, which coincides with our result, but comes with 
a larger uncertainty. As for the isospin asymmetry in the K* system, Eq. (70), we obtain 
a slightly lower value than Kagan and Neubert, Ref. [20], which is mainly due to their 
lower value A^^ = 0.35 GeV, compared to 0.51 GeV used by us, see Tab. 3. In the pQCD 
approach, the quoted asymmetry is about half of ours, but comes with a similar relative 
uncertainty [15], so that we agree within errors. Concerning, finally, the time- dependent 
CP-asymmetry S{K*^) in (76), we find approximate numerical agreement with Ref. [15], 
where the quark loops were modeled by intermediate vector states. As emphasized earlier, 
the exact size of the quark-loop contributions in this channel is not crucial since it is small 
compared to the leading term in ms/m;, sin (2/3), Eq. (77). The CP asymmetries S{pj) and 
S'(a;7) were also calculated by AP, but unfortunately their formulae miss the very crucial 
point that in S — > one has to deal with two physically distinguishable final states, 
namely Vl^r and Vr7l, whose amplitudes must be added incoherently, not coherently as 
done in Refs. [9, 10]. We therefore refrain from a direct comparison with their results. 

As for the relevance of our results for NP searches, the time-dependent CP asymmetries 
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are the cleanest observables since they are very small in the SM and constitute "quasi null 
tests" of the SM [67], in the sense that any measurement of a significantly non-zero value of 
these observables will be an unambiguous signal of NP. For K* , the asymmetry has already 
been measured, but is compatible with zero within errors. The asymmetry in Bg — > 07 is a 
very promising observable for the LHCb. Also the isospin asymmetry Ai{K*), Eq. (70), is 
very interesting for NP searches, and would become even more interesting upon completion 
of the NLO calculation started by Kagan and Neubert, Ref. [20], by including in particular 
the radiative corrections to the annihilation contribution shown in Fig. 1, with the photon 
emitted from the final-state quark lines. In contrast, neither the isospin asymmetry 
between p° and nor the asymmetry between and u are likely to be sensitive to 
NP. As for the branching ratios, we have, motivated by the inclusive B result, 
assumed no significant NP effects m B ^ K*^, and as long as there is no breakthrough 
in the calculation of the absolute values of form factors, any moderate NP effects in the 
branching ratios are likely to be obscured by the uncertainties. 

In summary we feel that exclusive b {s,d)^ transitions have a massive discovery 
potential for NP and envisage a great future at the LHC, which may be surpassed only 
by that of 6 ^ (s, o?) decays. 

Note added. After completion of the calculations presented in this paper, Ref. [44] 
appeared which contains a lattice calculation of Tf^^*{0) and in the quenched ap- 
proximation. The results are ^-^'(O) = 0.24 ± 0.031°;°^ and = 1.2 ± 0.1. The latter 
agrees with ours, Eq. (55), but comes with a slightly larger central value and uncertainty, 
while the former is a bit on the low side of the LCSR prediction given in Tab. 3. 
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A Vector Meson Decay Constants Revisited 

There are two types of decay constants for vector mesons: the vector coupling /y, for a 
longitudinally polarised meson, and the tensor coupling fy, for a transversely polarised 
meson: 

{Q\qii,D\V(p, e)) = ep^mvfv , {0\q(^nuD\V(p, e))^ = i(e^p^ - e^p^)fy{ij) . (A.l) 

Note that fvif^) depends on the renormalisation scale. 

The numerical values of these couplings are essential for our calculations. Whereas 
the extraction of the charged mesons' vector couplings from experimental data is straight- 
forward, that of the neutral mesons' p°, cu and is complicated by the mixing of these 
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particles and deserves a more detailed discussion, which we will give in Sec. A.l. The 
tensor couplings are not accessible experimentally, but have to be determined by non- 
perturbative methods, for instance QCD sum rules and lattice simulations. In Sees. A. 2 
and A.3, we briefly review the most recent results from these calculations. 



A.l Longitudinal Decay Constants from Experiment 

A. 1.1 The Charged Decay Constants fp-,K*- from r Decays 

The longitudinal decay constants of charged vector mesons can be extracted from t' 
V''i'r, with the measured branching ratios [28]: 

B{t- p-Ur) = (25.2 ± 0.4) X 10"^ , B{t~ K*-Ur) = (1-29 ± 0.05) x lO'^ . 

The decay rate is given by 

r(r- - V-u^) = iK^r/^ fl - ( 1 + 2^"^ 



IGtT -f I "-^ I \^ ^^2 J y ^2 

With \Vud\ = 0.9738 ± 0.0002 and |\4s| = 0-227 ± 0.001 [28], we get 

fp- = (210 ± 2b ± Ir.) MeV , Jk*- = (220 ± 4^ ± Ir. ± l|v-„^|) MeV , (A.2) 

where we have taken into account the uncertainties in the branching ratios, total decay 
rates and CKM matrix elements. The uncertainties of other input parameters are irrel- 
evant and the size of neglected corrections to the decay rate in a and higher powers in 
l/rri^r is expected to be smaller than the total uncertainty. 

A. 1.2 The Neutral Decay Constants fpo,u},4> from V" e+e~ 

The decay constants of p", uj and can be extracted from the electromagnetic annihilation 
process e+e^, which is, however, complicated by the mixing of these mesons. The 

states of definite isospin are given by 

= -^(|m) - \dd)) , la;,) = ^{\uu) + \dd)) , 10,) = \ss) , (A.3) 

where p has isospin 1 and uj and have isospin 0. In view of the subtleties of mixing, let us 
state clearly that the neutral decay constants /^o, and shall denote the coupling of 
the real particles to their isospin currents, e.g. {Q\s^^s\(l)) = e^m^f^. p-oj mixing violates 
isospin and hence is a purely electromagnetic effect which can be parametrized as 

\p) ~ \pi) - ep^\ui) , \u) ~ + ep^\pi) (A.4) 

with = 5p^/{{m^ - iT^/2f - (m^ - %Tp/2f) and 5^^ = -(0.004 ± 0.002) GcV' [69], 
which results in = (0.036±0.018)ie°'^^*; since rUp ~ m(^, e^^^ is almost purely imaginary. 
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This parameter was also determined experimentally [70] . The mixing of u and 0, on the 
other hand, is due to strong interactions: 

\uj) ~ \uji) - €^^\(f)i) , \(f)) ~ \(f)i) + e^^lcji) . (A.5) 

The mixing parameter has been determined to be e^^^ = 0.045 ±0.01 [71] by parametrising 
the SU(3) breaking in order to match the Gell-Mann-Okubo mass relation for light mesons 
with the observed masses. A direct measurement of this quantity was reported in a; — > 
e+e~ decays [72]. Evidently a full description of the mixing would involve all three states, 
but p-(f) mixing is expected to be very small because it is a second-order effect that requires 
both electromagnetic and strong interactions to be at work. 

As the — > e'^e" transition is an electromagnetic decay process, one needs the 
relevant electromagnetic currents of light quarks: 

jT = QuUl^^u + Qjj^d + g,s7^s = ^ (j;=° + 3 j'=') - hj^s ; (A.6) 

the isospin currents are defined as j^T^^^ — ■^(u'fuU ± d'-fij,d). The experimental rates are 
[28]: 

^ e+e-) = (4.7 ± 0.08) x 10"^ , B(u; e+e~) = (7.18 ± 0.12) x 10"^ , 
S(0 ^ e+e") = (2.97 ± 0.04) x 10"^ . (A.7) 
The theoretical expression for the decay rate is given by 

r{V'^e^e-)^^^fyCy, (A.8) 

3 my 

where the coefficients cy in the limit of no mixing can be read-off from (A. 3) and (A.6): 
Cpo = {Qu - Qdfll = 1/2, c^, = {Qu + Qdy/2 = 1/18 and c^, = = 1/9. The effect 
of, for instance, p-u mixing is to change CpO to CpO = \y^c~a — ep^^^fc^j\^ l\\ + epj(^ and 
correspondingly for c^. Including the mixing effects we finally get 

/po = (222±2Br±lrjMeV, 

= (187 ± 2Br ± lr„ ± 4,,^ ± MeV , 

= (215 ± 2Br ± Ir, ± 4^<^) MeV , (A.9) 

where again the uncertainties in the other input parameters are irrelevant and the cor- 
rections to (A.8) are expected to be smaller than the total uncertainty, p-io mixing has a 
negligible effect on /pO , but raises /^j by 2 to 3 MeV. a;-0 mixing is much more relevant 
and lowers by about 10 MeV and by about 13 MeV . 
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A. 2 Decay Constants from QCD Sum Rules 



The calculation of decay constants from QCD sum rules was one of the earhest appli- 
cations of this method [73]. More recent determinations include more (radiative and 
mass) corrections and updated values of input parameters. The most recent results were 
obtained in Refs. [22, 34] and read 

fp =(206±7)MeV, f^{lGeV) =(165±9)MeV, 

/i^* = (222 ± 8) MeV , /j^. (1 GeV) = (185 ± 10) MeV . (A.IO) 

Note that the determination of is more complicated than that of the other couplings 
and requires the inchision of higher resonances in the hadronic dispersion relation [34]. 
The above results refer to charged mesons; isospin breaking and meson mixing are not 
included. For comparison with lattice results, it proves convenient to also quote the results 
for the ratio of couplings: 

^1 (2 GeV) =0.70 ±0.04, [^J (2 GeV) = 0.73 ± 0.04 . (A.ll) 
A. 3 Decay Constants from Lattice QCD 

The ratio of decay constants fy/fv has been calculated by two lattice collaborations, in 
the quenched approximation. Ref. [74] obtains 

(2 GeV) = 0.72 ± 0.02 , ( ^ ) (2 GeV) = 0.74 ± 0.02 , 

\fK*J latt 

(2 GeV) = 0.76 ± 0.01 , (A.12) 
in the continuum limit, whereas Ref. [75] quotes: 

^1 (2 GeV) = 0.742 ± 0.014 , (^| (2 GeV) = 0.780 ± 0.008 , (A.13) 

Jp / laXt \f(t> J latt 

at the finite lattice spacing a = O.lOfm. 

The results from both lattice collaborations are roughly in agreement. It is evident 
that the ratios depend only weakly on the quark masses. 





1 


nj 


latt 




1 




latt 



A. 4 Discussion, Conclusions and Results 

In this paper we use the experimental results (A. 2) and (A. 9) for the longitudinal decay 
constants, averaging the two results for the p meson. For the tensor couplings of p and 

i°Including the NLO scaling factor /^(2GeV)//^(l GeV) = 0.876. 
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K* we use the sum rule results (A. 10). For a;, we assume isospin symmetry of the ratio 
of decay constants and use /J- (2 GeV) //^^ = 0.71 ± 0.03, which is the average of QCD 
sum rule and lattice results, to obtain a value for /J" from the measured f^. Finally, for 
(f) we use the lattice ratio (A. 12), with the more conservative uncertainty ±0.03, and the 
experimental value for Our final results which enter Tab. 3 are: 



The experimental vector couplings, (A. 2) and (A. 9), come with rather small uncer- 
tainties and indicate an isospin breaking in /pO,± of ~ 5%. Is it really justified to average 
both results into only one decay constant for the p? In order to answer this question, 
let us have a look at the experimental information on isospin breaking in other light- 
meson decay constants. For the tt, PDG gives /yri = (130.7 ± 0.4) MeV, whereas /^o is 
extracted from 7r° — > 77 as (130 ± 5) MeV which is perfectly compatible with /^±; the 
uncertainty is dominated by that of the 7r° lifetime [76]. A further confirmation of the 
smallness of isospin breaking comes from the CLEO measurements of — > {K~ , T^~)e^v 
and — > (K°, 7r°)e''"i/ [45]. These decays are sensitive to the tt and K meson decay 
constants via the form factors, which, at least in the LCSR approach, are directly pro- 
portional to f-K,K-, see e.g. Ref. [77]. The data for D Kev indicate that the isospin 
breaking of the form factor is (3 ± 2)%. Although this result also includes potential 
isospin breaking of both ]d and the dynamical part of the form factor, the corresponding 
effects are neither expected to be sizable nor to cancel each other, so that indeed for both 
TT and K decay constants isospin breaking is smaller than 5% at \a. Taking this as an 
indication for the generic size of isospin breaking in light mesons, we conclude that a 5% 
difference between /^o and jp± is not excluded, but on the large side. This implies that 
it is indeed appropriate to average the experimental results for neutral and charged p as 
done in (A. 14). 

Note that the QCD sum rule results for the vector couplings, (A. 10), agree rather 
well with the experimental results, (A. 2) and (A. 9), which increases confidence in the 
corresponding results for the tensor couplings, particularly as the sum rule results for the 
ratios, (A. 11), also agree with those from lattice, (A. 12). Nevertheless there is one fact 
which remains somewhat puzzling, namely the difference of nearly 20% between the lo 
and (P couplings from — > e^e~ ^ which is larger than the expectation from QCD sum 
rules. This difference could be caused by electromagnetic corrections, different values 
of the up and down quark condensate and differences in the continuum thresholds and 
Borcl windows. The latter effects should not exceed the typical accuracy of sum rules 
themselves, which is about 10%, and the former effects are expected to be very small. 
On the other hand, the individual experimental results entering the PDG averages (A. 7) 




(216 ± 3) MeV , 
(187 ±5) MeV, 
(220 ± 5) MeV , 
(215 ±5) MeV, 



/p^(lGeV) 
/a^(lGeV) 
/i*(lGeV) 
/i(lGeV) 



(165 ± 9) MeV , 
(151 ± 9) MeV , 
(185 ±10) MeV, 
(186 ±9) MeV. 



(A.14) 
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are spread over a wide range, particularly for a;, which indicates that the uncertainties 
quoted in (A. 7) may be on the optimistic side. 

There remains one subtle point to be discussed, namely that for LCSRs for b ^ dj 
transitions and p° or u in the final state, one needs the decay constant 

V2{0\d^^d\V''ie)) = e^mvo f'^} , 

rather than /yo. The quantity fyo could differ from /yo through mixing with the other 
neutral mesons. Fortunately, uj-4> mixing is irrelevant because the coupling of (j) to the 
down-quark current is highly suppressed, and p-uj mixing has a small effect because the 
mixing parameter is almost imaginary. The total impact of mixing is hence below the 
theoretical uncertainty and can safely be neglected. 
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